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I.   INTRODUCTION 

A  Foundation  Deformation  Prediction  Symposium  was  held 
at  MIT  during  November  13  -  15,  1974.   The  purpose  of  the 
symposium  was  to  examine  current  design  methodology  through 
the  comparison  of  predicted  and  measured  performance  of  a 
soft  soil  foundation  subjected  to  an  additional  embankment 
load. 

Since  1967  the  Department  of  Civil  Engineering,  Massa- 
chusetts Institute  of  Technology,  has  participated  in  research 
on  the  performance  of  a  high  embankment  placed  over  a  thick 
clay  foundation.   The  Massachusetts  Department  of  Public  Works 
and  the  Federal  Highway  Administration  have  sponsored  the  re- 
search.  The  major  objectives  of  the  research  have  been: 

1.  To  develop  methods  for  predicting  the  performance 
of  the  clay  foundation  under  an  embankment  load; 

2.  To  measure  the  performance  of  the  foundation 
during  and  after  construction,  and 

3.  To  compare  the  measured  and  predicted  performances. 

A  substantial  amount  of  field  performance  information  has 
been  obtained  and  analyzed. 

As  a  part  of  this  research,  additional  fill  was  placed  on 
the  present  embankment  during  late  summer  19  74  to  cause  large 
undrained  deformations  in  the  clay  foundation  (i.e.,  reaching 
limiting  equilibrium)  and  measurements  of  deformations  and  pore 
pressures  were  taken.   Ten  engineers  (or  teams)  submitted 


predictions  of  key  aspects  of  performance  to  two  referees 
prior  to  the  conduct  of  the  field  test. 

All  of  the  predictors  were  supplied  the  same  basic  infor- 
mation.  Several  predictors  requested  additional  data. 
Whenever  possible,  these  additional  data  were  then  supplied 
to  all  predictors.   All  predictions  were  given  to  an  inde- 
pendent referee  before  the  start  of  the  field  test.   No  pre- 
dicted or  measured  values  were  revealed  until  the  symposium. 

The  symposium  brought  together  the  group  of  predictors  to  focus  their 
attention  on  a  single  problem  rather  than  present  an  analysis  of  different 
problems.  In  this  way,  the  maximum  synergistic  benefit  could  be  obtained 

from  the  efforts  of  the  group. 

At  the  symposium,  each  predictor  presented  his  approach 

to  predicting  performance  at  previously  specified  locations 

and  his  predictions.   Measured  values  were  then  revealed. 

Each  predictor  commented  on  his  predictions  and  approach  in 

light  of  the  field  performance.   Discussion  of  the  predictions 

and  measured  performance  by  the  audience  followed. 

These  proceedings  document  the  symposium  including  pre- 
dictions, measurements  and  discussion.   Written  presentations 
submitted  by  the  predictors  are  presented  where  available; 
otherwise,  we  have  attempted  to  indicate  the  approach  taken, 
the  parameters  used  and  the  predicted  values  as  recorded  dur- 
ing the  symposium.   In  the  interest  of  preserving  the  flavor 
and  character  of  the  symposium,  we  have  not  tried  to  compile 
a  polished  set  of  proceedings. 

Prof.  T.  W.  Lambe  opened  the  symposium  by  discussing  the 
importance  of  predictions.   Geotechnical  work  at  MIT  during 
the  last  two  decades  has  led  to  formulation  of  the  following 
important  considerations: 


1.  Decisions  are  based  on  predictions. 

2.  Predictions  are  an  important  part  of  civil 
engineering. 

3.  Predictions  must  be  evaluated  to  assist  in  the 
decision  making  process  and  to  improve  the  pre- 
diction process. 

4.  Predictions  are  frequently  inaccurate. 

5.  The  major  cause  of  inaccurate  predictions  is 
faulty  and  insufficient  data. 

Several  key  steps  can  be  identified  in  the  prediction 
process: 

1.  Determine  the  field  situation  (gather  facts); 

2.  Simplify  the  field  situation  (organize  facts); 

3.  Determine  mechanism (s); 

4.  Select  both  Method (s)  and  Parameters  ; 

5.  Manipulate  Method  and  Parameters — make  predic- 
tion; 

6 .  Portray  prediction,  and 

7.  Evaluate- Improve  (during  and  after  symposium). 

To  help  classify  predictions,  we  at  MIT  use  the  following 
descriptions: 

Prediction  When  prediction  Results  at  time 

type  made  prediction  made 

Type  A  Before  event  

Type  B  During  event  Not  known 

Type  C  After  event  Not  known 
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This  symposium  was  organized  to  provide  an  opportunity 
to  present,  discuss  and  evaluate  Type  A  predictions. 


NOTE 


The  original  instructions  to  each  of  the  predictors 
did  not  call  for  the  preparation  and  submission  of  a 
written  report  describing  his  prediction.   This  fact 
explains  the  variation  in  the  detail  and  form  of  the 
contributions.   For  the  sake  of  accuracy  we  chose  not 
to  redraft  certain  of  the  figures. 
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II.   PROBLEM  STATEMENT 

Embankment  Description 

Figure  II-l  shows  a  cross-section  of  the  embankment  con- 
structed in  1967  to  1969  and  the  subsurface  soil  profile  at 
Station  263.   Plans  were  to  widen  a  300-foot  long  section  of 
the  embankment  crest  centered  at  Station  263  by  dumping  fill 
along  the  east  slope   until  the  slope  was  about  40°.   Trucks 
were  then  to  dump  fill  on  top  of  the  embankment  in  1-foot  lifts 
at  an  average  rate  of  1-  to  2-foot  height  increase  per  day, 
excepting  weekends.   No  compaction,  other  than  that  produced 
by  normal  construction  traffic,  was  planned.   The  additional 
fill  would  come  from  a  portion  of  the  present  embankment  about 
1800  feet  away.   Filling  was  to  stop  for  a  two-day  period  after 
6  feet  of  additional  fill  had  been  placed  on  the  embankment. 

Predictions  to  be  Made 

The  predictors  were  asked  to  predict  the  following  quan- 
tities : 

1.  After  6  feet  of  fill  have  been  added  to  1969  embank- 
ment, predict: 

a.  additional  horizontal  movement  of  SI-3  and 
SI-4  at  elevations   -30  and  -70; 

b.  additional  settlement  of  SP-1; 

c.  additional  pore  pressure  at  P-3,  P-4,  P-6 
(maximum  during  two-day  period  after  filling 
is  stopped) t    and 
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d.   heave  at  H-l,  H-2. 

2.  Additional  height  of  fill  (T  in  Fig.  II-l) 
necessary  to  cause  a  stability  failure  (dif- 
ferential deformation  at  crest  >  0 . 5  m  con- 
stitutes failure) . 

3.  Pore  pressure  at  P-3,  P-4,  P-6  at  failure. 
Figure  II-l  shows  instrument  locations. 
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III.   DATA  SUPPLIED  TO  PREDICTORS 

Information  and  data  were  sent  to  each  predictor  in  a 
series  of  bulletins.   This  section  summarizes  the  information 
given  in  these  bulletins. 

History  of  Initial  Embankment  Construction 

Initial  construction  began   in  August,  196  7  with  removal 
of  the  top  layer  of  silt  and  peat  and  placement  of  a  10-foot- 
thick  sand  and  gravel  working  mat.   During  placement  the  fill 
was  end-dumped  with  no  compaction  up  to  elevation  +5  ft.   Above 
elevation  +5  ft,  the  fill  was  placed  uniformly  along  the  entire 
length  of  the  embankment  with  scrapers  and  trucks.   Rubber 
tired  tractors  and  two  8,000  gallon  water  tanks  (filled  with  salt 
water)  helped  compact  the  fill  to  an  average  total  unit  weight 
of  119   lb/ft3. 

The  sequence  of  construction  was: 

August  196  7  —  Peat  removal  and  replacement  commended. 

April   19  6  8  —  Peat  removal  and  replacement  completed. 
Embankment  filling  commenced. 

November  1968-- All  filling  stopped  (Nov.  15). 

April   1969  — Filling  resumed  (April.  15). 

May     1969  — Filling  100%  complete. 

Geology  of  Area 

The  5-  to  45-foot-thick  layer  of  glacial  till  found  at 
the  site  was  probably  deposited  during  the  Wisconsin  glaciation 
over  a  gray  shale  bedrock  found  at  depths  of  70  to  204  feet. 
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Radiocarbon  dates  indicate  that  the  Boston  area  became  free  of 
glacial  ice  at  least  14,000  years  ago  and  deposition  of  the 
clay  was  in  progress  at  that  time.   Sediments  of  glacial  origin 
were  deposited  in  brackish  water.   Subsequently,  during  the 
Valders  glacial  substage  (12,250  to  11,750  years  ago),  sea  level 
fell  with  respect  to  land.   The  clay  sediments  emerged  from  be- 
low sea  level  and  were  eroded  by  streams  and  weathered  sub- 
aerially.   The  weathering  and  subsequent  desiccation  resulted 
in  the  formation  of  a  stiff  "crust"  which  exhibits  the  properties 
of  overconsolidation.   The  Valder's  glacial  substage  was  followed 
by  warmer  climates  during  which  sea  level  rose  rapidly  and  sands 
were  deposited  in  depressions  in  the  clay  surface.   As  the  sea 
continued  to  rise  at  a  slower   rate  to  its  present  level,  the 
entire  area  was  covered  by  organic  silt,  shells  and  peat. 

Preconstruction  Soil  Properties 

The  project  lies  in  a  low  tidal  marsh  which  is  subjected 
to  daily  tidal  inundations.   The  average  ground  elevation  be- 
fore construction  was  +5  feet  (USGS) .   The  average  daily  tide 
variation  is  2.5  feet  with  average  tide   elevation  at  +2.5  ft 
(USCS) . 

Preconstruction  soil  properties  of  the  clay  from  labora- 
tory and  field  tests  at  nearby  Station  246  are  shown  in  Figure 
III-l.   Very  little  preconstruction  testing  was  done  at  Station 
263. 


The  thickness  of  the  natural  sand  layer  overlying  the 
clay  varies  with  location.   On  the  west  side  of  Station  263, 
the  sand  grades  into  clay  at  elevation  -12  feet.   Borings  at 
the  centerline  of  Station  263  indicate  the  sand-clay  interface 
is  presently  (1973)  at  elev.  -18  to  -20  feet  while  on  the  east - 
side  sand  ends  at  elevation  -7  feet.   Blow  counts  from  standard 
penetration  tests  in  the  sand  at  Station  263  vary  from  10  to 
49.   At  Station  246  the  sand  blow  counts  are  8  to  25  with  an 
average  of  17.   No  laboratory  tests  were  run  on  this  material. 
Split  spoon  samples  indicate  that  it  is  a  well-graded,  silty  sand 
with  some  gravel. 

No  laboratory  tests  were  run  on  the  underlying  till.   It 
is  quite  dense  with  blow  counts  generally  in  excess  of  50  but 
varying  from  24  to  171.   The  following  elevations  have  been 
determined  for  the  top  of  the  till  near  Station  263: 


Location 

Elevation 

(ft) 

SI-1 

-109 

SI-2 

-108 

SI-3 

-87 

SI-4 

-65 

140  feet  east  of  g 

-68 

At  toe  Sta.  263  +90 

-80 

At  toe  Sta.  263  +10 

-87 
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The  fill  consists  of  well-graded  fine  to  coarse  sand, 
with  some  fine  to  medium  gravel  and  a  few  boulders.   The  pre- 
dominant mineral  is  angular  quartz,  but  some  mica  is  present. 
Figure  III-2  shows  the  results  of  three  drained,  triaxial 
loading  tests  on  the  fill  material. 

1974   Soil  Properties 

A  series  of  laboratory  and  field  tests  were  conducted  on 
the  clay  at  Station  263  just  prior  to  the  addition  of  fill  for 
the  field  test.  Figures  III-3  and  III-4  show  a  summary  of  the 
results  of  these  tests.  All  lab  tests  were  run  on  samples  trimmed 
from  5-inch  diameter  undisturbed  samples  obtained  a  few  months 
prior  to  the  start  of  the  field  test. 

Field  Measurements  of  Previous  Performance 

Three  piezometers,  four  slope  indicators  and  one  settle- 
ment plate  were  installed  during  the  early  stages  of  construction. 
The  three  piezometers  at  the  centerline  (A,  B,  and  C)  ceased  to 
function  by  1971  and  were  replaced  by  piezometers  P3,  P4,  and 
P5  in  July,  1973,  along  with  additional  piezometers  to  monitor 
embankment  performance  during  additional  fill  placement.   The 
locations  of  these  instruments  are  given  in  Figure  II-l  and  the 
measured  results  are  sumarized  in  Figure  III-5. 

Publications  Describing  1-95  Research 

The  following  is  a  list  of  professional  papers  which  dis- 
cuss aspects  of  1-95: 
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D'Appolonia,    D.J.,  T.W.  Lambe  and  H.G.  Poulos    (1971)    "Evaluation  of 
pPo?e   Pressures   Beneath    an    Embankment",    JSMFD  ASCE,    Vol.    97 


SM    6. 


Lambe,  T.W.  (1973)  "Predictions  in  Soil  Engineering",  13th 
Rankine  Lecture,  Geotechnigue,  Vol.  23,  No.  2. 

Lambe,  T.W. ,  D.J.  D'Appolonia,  K.  Karlsrud  and  R.C.  Kirby 
(1972)  "The  Performance  of  the  Foundation  Under  a  Hign 
Embankment",  Journal  of  the  Boston  Society  of  Civil 
Engineers,  Vol.  59,  No.  2,  April. 


Additionally,  there  are  several  MIT  research  reports  document- 
ing various  aspects  of  the  joint  research  effort  and  a  series 
of  bulletins  giving  more  information  on  the  soil  properties  at 
Sta.  263. 
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IV.   PRESENTATION  OF  PREDICTIONS 

General 

Field  loading  took  place  after  the  formal  predictions 
were  submitted,  but  prior  to  the  symposium.   Before  disclosing 
the  field  results,  each  presented  his  predictions  along  with 
a  discussion  of  the  methods  he  employed.   Several  of  the  pre- 
dictors used  normalized  soil  properties  (Ladd  &  Foott,  1974) 
and  therefore,  a  discussion  of  this  concept  was  presented  by 
Prof.  C.  C.  Ladd  (one  of  the  predictors)  prior  to  the  presen- 
tation of  the  predictions.   A  short  biographical  note  for  each 
predictor  is  given  in  Section  VII. 

Several  predictors  provided  papers  documenting  their  pre- 
dictions.  These  papers  are  in  Appendices  A  through  H.   The 
following  pac,_  present  summaries  of  all  the  predictions. 

H.  G.  Poulos  and  E.  H.  Davis   (See  Appendix  A  also) 

Davis  and  Poulos  based  their  prediction  on  two  finite 
element  computer  analyses: 

(1)  An  elasto-plastic  analysis  was  employed  to  deter- 
mine the  undrained  behavior  of  the  embankment 
during  construction.   The  finite  element  mesh  for 
this  analysis  is  shown  in  Fig.  IV- 1. 

(2)  An  elastic  consolidation  analysis  (Biot  Theory) 
was  used  to  estimate  pore  pressures  and  consoli- 
dation settlements  both  during  and  after  construc- 
tion. 
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In  addition,  several  independent  approximations  and  hand 
calculations  were  carried  out  to  check  the  finite  element  cal- 
culations. 

The  prediction  process  was  divided  into  two  stages. 
Stage  I  predicted  the  performance  of  the  embankment  to  its  1969 
height  for  which  field  measurements  were  available.   Stage  II 
predicted  the  behavior  of  the  embankment  resulting  from  the 
recent  addition  of  fill. 

The  soil  parameters  for  the  Stage  I  finite  element  analy- 
ses were  originally  adopted  from  laboratory  and  field  test  re- 
sults  and  published  data.   Some  of  them  were  adjusted,  however, 
by  trial  and  error  during  Stage  I  so  that  the  resulting  predict- 
ions agreed  reasonably  well  with  actual  observations.   The  range 
of  values  of  all  soil  parameters  are  shown  in  Table  IV-1. 

In  Stage  II,  allowance  was  made  for  the  effects  on  consoli- 
dation which  had  occurred  after  1969,  but  prior  to  the  1974 
addition  of  fill,  by  determining,  from  the  results  of  the  finite 
element  consolidation  analysis,  the  changes  in  effective  stress 
which  had  occurred  since  the  start  of  embankment  construction. 
The  values  of  s   and  E   just  prior  to  the  addition  of  fill  were 
estimated  on  the  basis  of  the  current  vertical  effective  stress, 
and  hence  the  current  OCR,  and  the  correlations  with  OCR  of  the 

ratios  s  /c    and  E   /a    established  in  Stage  I. 

u  '  vo       u  /  vo  3 

The  contours  of  the  estimated  values  of  s  and  E  before 
the  19  74  addition  of  fill,  as  shown  in  Figures  IV-2  and  IV- 3, 
were  used  in  the  finite  element,  elasto-plastic,  undrained 
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analysis  for  the  case  of  6-feet  additional  fill.   Predictions 
of  deformations  and  additional  pore  pressures  (calculated  as 
the  mean  of  the  vertical  and  the  horizontal  stress  increases) 
were  obtained. 

Further  element  layers  were  then  added  to  the  top  of  the 
embankment  and  general  failure  was  assumed  to  occur  when  the 
rate  of  increase  of  plastic  work  became  negative  and  convergence 
could  not  be  obtained.   The  resulting  predictions  of  the  addi- 
tional height  of  fill  to  cause  stability  failure  and  the  other 
required  predictions  are  listed  in  Table  IV-2. 

The  UU  test  results  were  disregarded  and  the  field  vane 
test  results  were  not  directly  used.     Strength  anisotropy  of 
the  Boston  Blue  Clay  was  considered  significant  and  therefore, 
the  average  s  was  taken  to  be  90  percent  of  the  average  of  s 
in  the  vertical  and  horizontal  directions.   The  Stage  I  analysis 
also  led  the  predictors  to  conclude  that:  (1)  undrained  creep 
could  be  a  major  problem  in  the  Boston  Blue  Clay,  (2)  the  till 
beneath  the  clay  was  effectively  impermeable,  and  (3)  hand  cal- 
culation methods  of  approximation  could  give  very  good  results, 
if  accurate  values  of  soil  parameters  were  used. 

W.  W.  Moore 

Serving  as  head  of  the  prediction  team  for  Dames  &  Moore, 
Mr.  Moore  opened  the  presentation  with  some  comments  on  the  role 
of  geotechnical  engineering  consultation  services  and  his  group's 
eagerness  to  participate  in  this  event.   Mr.  M.  J.  Hess  then 
took  over,  introduced  the  rest  of  the  team~-Mr.  R.  A.  Bell  and 
Mr.  W.  J.  Roberds — and  presented  the  analytical  procedures 
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involved  and  the  resulting  Dames  &  Moore's  predictions.   The 
predictions  developed  by  Dames  &  Moore  used  all  of  the  subsurface 
information  that  was  provided,  as  well  as  the  various  published 
articles  on  research  work  performed  by  MIT  on  the  Interstate 
9  5  embankment. 

Strength   properties  of  the  Boston  Blue  Clay  were  assessed 
by  the  SHANSEP  approach  developed  by  Professor  C.  C.  Ladd  of 
MIT.   Stresses  induced  by  the  embankment  were  computed  by  elas- 
tic theory.   Pore  pressure  distribution  was  estimated  assuming 
that  change  in  pore  pressure  was  equal  to  the  change  in  octa- 
hedral normal  stress,  up  to  the  point  of  local  yield  and  at 
higher  stresses  the  change  in  pore  pressures  was  assumed  equal 
to  the  major  principal  stress  increment.   The  estimated  values 
were  compared  with  previous  piezometer  readings  at  the  test  sec- 
tion.  The  dissipation  of  pore  pressure  was  assumed  by  using  a 

2  2 

coefficient  of  consolidation  of  1.44  feet   per  day  and  0.14  ft 

per  day  for  overconsolidated  and  normally  consolidated  Boston 
Blue  Clay,  respectively  for  the  stress  ranges  concerned  in  these 
analyses.   The  higher  value  for  the  overconsolidated  clay  applies 
only  to  stresses  below  the  preconsolidation  pressure.   Drainage 
was  assumed  to  occur  at  both  the  top  and  bottom  of  the  layer. 
A  profile  of  excess  pore  pressure  existing  as  of  August  1974 
was  developed  from  piezometer  readings  taken  in  the  field  at 
that  time,  as  shown  in  the  attached  Figure  IV-4.   With  these 
values  of  excess  pore  pressures,  and  the  elastic  stress  analysis 
due  to  the  embankment  loading,  a  profile  of  vertical  effective 
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stresses  was  developed,  as  shown  in  Figure  IV-5.   By  using  the 
effective  stress  distribution   and  the  estimated  maximum  past 
vertical  effective  stresses  (Figure  IV-6)  ,  the  shear  strength  pro- 
file as  shown  in  Figure  IV  7,  was  plotted  by  using  SHANSEP  theory.  The 
strength  of  the  fill  was  evaluated  assuming  cf>  =  40°  in  the  existing  embank- 
ment and  <$>  =35°  in  newly  placed  fill.  The  normal  stress  on  the  portion  of 
failure  surface  in  the  fill  was  computed  by  elastic  theory  and  also 
considered  that  an  active  state  of  stress  existed.   Stability 
analyses  by  both  circular  arc  and  sliding  wedge  were  performed 
with  the  previously  described  strength  parameters.   Figure  IV-7 
shows  the  critical  surfaces  for  both  modes  of  failure.   The 
circular  arc  analyses  resulted  in  a  factor  of  safety  equal  to 
1.0  with  the  top  of  fill  at  elevation  +55  ft.   The  results  of 
the  circular  arc  analysis  were  largely  influenced  by  the  strength 
of  the  fill.   With  the  sliding  wedge  analysis,  a  factor  of  safety 
of  1.0  occurred  with  the  top  of  fill  at  elevation  +53  ft. 

The  creep  behavior  of  the  Boston  Blue  Clay  was  assessed  in 
a  qualitative  manner  when  reviewing  the  results  of  the  stability 
analysis.   The  assessment  was  based  on  past  experiences  of  slope 
failures  in  San  Francisco  Bay  mud.   Even  though  the  loading  of 
the  test  embankment  was  to  be  in  a  relavitely  short  period,  the 
effects  of  a  possible  delay  were  considered  in  an  alternate  pre- 
diction. 

Failure  of  the  embankment,  if  continuously  loaded  at  a  rate 
of  1.5  ft  per  day,  was  predicted  at  elevation  +53.   If  loading 
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was  delayed  on  the  order  of  5  to  7  days,  failure  was  predicted 
at  elevation  +49  to  +50.   The  latter  prediction  estimated  the 
available  shear  strength  of  the  Boston  Blue  Clay  to  be  on  the 
order  of  80  percent  of  the  ultimate.   Laboratory  tests  to  sim- 
ulate the  varying  levels  of  field  stress  under  variable  dura- 
tion of  loading  were  considered  very  important  to  determining 
the  height  of  failure. 

Excess  pore  pressures  were  estimated  by  using  elastic 
theory,  with  total  applied  stresses  due  to  additional  fill. 
The  stress  distribution  charts  by  Giroud  &  Watissee  were  used 
with  a  Poissons  ratio  of  0.5.   Solutions  of  applied  total 
stresses  were  obtained  at  the  points  corresponding  to  piezo- 
meters P-3,  P-4,  and  P-6  for  the  two  loading  cases  of  an  addi- 
tional 6  ft  of  fill  and  for  an  estimated  15  feet  of  additional 

fill.   For  the  first  loading,  it  was  assumed  that  Au  =  Aa   .  , 

3  e    oct 

i.e.,  that  no  local  yielding  had  yet  occurred.   This  assumption 
was  found  to  be  consistent  with  the  Mohr's  circle  representa- 
tion of  the  stress  state.   The  additional  excess  pore  pressures 
for  an  additional  6  ft  of  fill  were  estimated  to  be  8.8,  7.2, 
and  4.8  feet  of  water  for  piezometers  P-3,  P-4,  and  P-6,  re- 
spectively.  Pore  pressures  for  an  additional  9  ft  of  fill 
(total  of  15  feet)  were  assumed  to  vary  with  Aa  .   At  the  pre- 
dicted failure  height,  these  values  were  estimated  to  be  16.6, 
16.6   and  11.9  feet  of  water. 


-24- 


Deformations  were  predicted  for  an  additional  6  feet  of 
fill,  which  was  kept  on  for  two  days.   The  additional  vertical 
settlement  at  point  SP-1  was  estimated  by  elastic  analysis, 
using  the  previously  developed  applied  stresses  and  strength 
profiles.   Assumptions  were  made  that  the  Boston  Blue  Clay  was 
linearly  elastic  within  the  loading  range  arid  that  E  =  1200 
S  .   The  additional  settlement  was  estimated  to  be  1  to  1  1/2 
inches. 

The  estimations  of  both  horizontal  displacement  and  heave 
were  based  primarily  on  the  slope  indicator  data  at  Station 
263  +  00  and  a  comparative  study  of  the  settlement,  heave,  and 
horizontal  displacement  measurements  taken  during  loading  of 
fuel  storage  tanks  at  a  soft  ground  site  in  Mississippi  and  from 
a  land  reclamation  project  at  the  Redwood  Shores  Development  in 
California.   In  the  comparative  studies,  rough  percentages  of 
heave  and  horizontal  displacement  were  compared  to  observed  set- 
tlements during  undrained  loading.   Horizontal  displacements 
were  thus  empirically  estimated  to  be  1/2  to  1  inch  at  elevation 
-30  and  less  than  1/2  inch  at  elevation  -  30  were  estimated  for 
SI-4.   Heave  at  15  to  30  feet  east  of  the  toe  of  the  slope  was 
estimated  to  be  less  than  1  inch  and  at  65  to  70  feet  east  of 
the  toe  to  be  less  than  1/2   an  inch. 

Table  IV-3  presents  the  predictions. 
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C.  C.  Ladd 

The  work  done  by  the  MIT  team  was  organized  and  directed 
by  Professor  C.  C.  Ladd.   MIT  had  access  to  considerable  back- 
ground data  for  the  1-95  site;  however,  their  prediction  was 
also  submitted  prior  to  the  beginning  of  loading.   Prof.  Ladd's 
predictions  are  summarized  in  the  following  paragraphs.   A  com- 
plete report  on  the  MIT  prediction  is  presented  in  Appendix  B. 

Prof.  Ladd's  predictions  were  based  on  Normalized  Soil 
Properties  for  Boston  Blue  Clay.   The  selection  of  strength  and 
deformation  parameters  for  the  analyses  is  treated  in  detail  in 
Appendix  B.   An  accurate  determination  of  the  clay's  maximum 
past  pressure  was  particularly  important  for  the  selection  of 
the  proper  normalized  parameters.   The  finite  element  program 
FEECON  was  utilized  for  the  deformation  and  pore  pressure  pre- 
dictions.  The  initial  construction  of  the  embankment  was  first 
analyzed  in  order  to  "calibrate"  the  computer  model  with  the 
measured  field  performance  presented  in  the  bulletins.   This  led 
to  relatively  minor  changes  in  values  of  the  moduli.   a  second  FEE- 
CON  analysis,  which  considered  changes  in  properties  due  to  par- 
tial consolidation  under  the  embankment,  was  performed  to  obtain 
the  required  predictions.   Figure  IV- 8  shows  the  input  methodology 
for  FEECON  2  and  defines  the  four  fill  stages  analyzed.   Figure 
IV- 9  shows  the  grid  in  the  area  of  interest  and  the  predicted 
zones  of  yielding  beneath  the  left  side  of  the  embankment  after 
adding  6  ft  of  fill  (end  of  Stage  III) . 

Figure  IV- 10  presents  computed  and  measured  deformation 
data  for  inclinometers  SI-3  and  4.   Line  II  is  the  measured 
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horizontal  movement,  ATT,  during  Stage  II  construction  (elev. 

n 

34  ft  to  40  ft).   Measured  increases  in  Au  over  a  2.5-year 

n 

period  at  constant  fill  height,  presumedly  due  to  creep,  are  denoted  by 
Line  C.  From  the  information  shown  in  Figure  IV-10,  it  was  concluded 
that  creep  deformations  were  very  significant  for  the  problem. 
Since  there  were  no  lateral  deformations  in  the  till,  the 
best  estimate  of  the  predicted  movements  was  taken  as  those 
computed  by  FEECON  2  at  Elev.  -30  ft  and  -37  ft  minus  those 
computed  in  the  till.   These  values  are  shown  at  the  bottom 
of  Figure  IV-10. 

Plots  of  FEECON  2  computed  heaves  at  the  ground  surface 
(top  of  peat  at  elev.  +5  ft)  and  at  the  top  of  the  clay  (elev. 
-10  ft)  vs.  distance  from  the  centerline  are  shown  in  Figure 
IV-11.   Two  elevations  were  considered  since  the  location  of 
the  bottom  of  the  heave  stakes  was  not  known.   The  best  esti- 
mate was  based  on  the  average  of  the  two  curves. 

Figure  IV-12  presents  a  plot  of  fill  elevation  vs.  cen- 
terline settlement  at  elev.  -31  ft  (original  top  of  natural 
sand)  and  elev.  -10  ft  (original  top  of  clay  computed  from  FEE- 
CON) .  The  predicted  values,  a  combination  of  initial  settlement 
due  to  undrained  shear  and  consolidation  settlement,  are  shown 
at  the  bottom  of  Figure  IV-12. 

Excess  pore  pressures  were  predicted  using  the  equation: 

Au  =   Aa  ,_      +     a   At 

oct  oct 

Table  IV-4  presents  the  predicted  values  of  Au  at  the  three  re- 
quired piezometer  locations  after  placing  6  feet  of  fill  on  the 
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existing  embankment.   A  range  of  values  of  "a"  were  used  based 
on  the  values  back-calculated  from  field  data  for  filling  from 
elev.  +35  ft  to  elev.  +40) ,   The  pore  pressures  at  failure 
were  computed  in  the  same  manner. 

Stability  analyses  performed  to  predict  the  added  height 
of  fill  to  cause  failure  concentrated  on  Morgenstern-Price 
(M-P)  analyses  using  normalized  values  of  undrained  strength 
s  .   Undrained  strength  anisotropy  was  taken  into  account  by 
varying  the  value  of  su  according  to  the  inclination  of  the 
failure  plane  as  shown  in  Figure  IV-13.   Some  circular  arc  an- 
alyses using  corrected  field  vane  strengths  were  made  as  a  semi- 
empirical  check  on  the  M-D  analyses.   The  results  from  the  sta- 
bility analyses  are  shown  in  Figure  IV-14. 

Professor  Ladd's  predictions  are  summarized  in  Table  IV-5. 


-28- 


Z.C.  Moh,  E.W.  Brand  and  A.M.  Richardson 

The  team  from  the  Asian  Institute  of  Technology  (AIT)  was 
represented  at  the  Symposium  by  Prof.  Z.  C.  Moh.   A  complete 
report  on  the  AIT  prediction  is  included  as  Appendix  C  .   Table 
IV- 6  lists  AIT's  predictions. 

The  additional  height  of  fill  necessary  to  cause  a  stab- 
ility failure  was  estimated  using  the  Swedish  slip  circle  method. 
Total  stress  analysis  was  used  since  it  was  assumed  that  the 
construction  would  be  carried  out  essentially  in  undrained  con- 
ditions. 

A  stability  analysis  was  performed  using  an  IBM  1130  digi- 
tal computer.   The  computer  programs  were  written  to  analyze 
the  stability  of  the  test  embankment  under  the  existing  condi- 
tions  taking  into  account  the  different  subsoil  strata  and  vari- 
ations of  strength  values  with  depth.   Variables  introduced  in 
the  analysis  were  mode  of  failure,  strength  values  of  the  foun- 
dation soils,  unit  weight  and  strength  parameters  of  the  embank- 
ment fill.   The  investigation  of  the  stability  of  the  existing 
embankment  considered  three  modes  of  failure: 

(1)  Assume  that  the  strength  of  the  embankment  and 
the  underlying  sand  layer  were  fully  mobilized 
along  the  circular  slip  surface; 

(2)  Assume  that  the  strengths  of  the  embankment  and 
the  underlying  sand  layer  along  the  slip  surface 
did  not  contribute  to  the  stability,  and 

(3)  Assume  the  strengths  of  the  embankment  and  the 
underlying  sand  layer  were  fully  mobilized  as 
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being  equal  to  the  active  earth  pressure  on 
a  vertical  plane  multiplied  by  tangent  of  c}> 
for  the  appropriate  material. 

The  analysis  for  conditions  at  failure  considered  two  modes  of 

failure: 

(1)  Assume  that  the  strengths  of  the  fill,  the 
embankment  and  the  underlying  sand  layer 
were  fully  mobilized  as  being  equal  to  the 
active  earth  pressure  on  a  vertical  plane 
multiplied  by  tangent  cf>  for  the  appropriate 
material  and 

(2)  Assume  that  the  strengths  of  the  fill,  the 
embankment  and  the  underlying  sand  were  not 
mobilized. 

Based  on  data  provided,  the  following  strength  parameters 
were  adopted: 

(1)  granular  embankment  fill:   <j>  =42°  from  drained  test  data. 

(2)  sand  layer  beneath  the  embankment:   <j)  =  38° 

Since  the  additional  fill  will  be  placed  in  a  loose  condition, 
the  effective  angle  of  shearing  resistance  was  taken  as  4>  =  36°. 
We  assume  that  the  sand  fill  will  have  a  placement  water  content 
of  10  to  12%  and  a  total  unit  weight  of  approximately  106  lbs 
per  cu  ft.   Bulletin  3  gave  the  total  unit  weight  of  the  exist- 
ing embankment  as  varying  from  102.1  to  134  lbs  per  cu  ft  with 
an  average  of  119  pcf. 

The  additional  fill  would  be  placed  in  a  short  period  of 
time.   Thus,  the  foundation  stability  analysis  was  considered 
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an  undrained  case.   The  total  stress  analysis  method  was  adopted. 
The  strength  parameters  chosen  were  from  field  vane  tests  and 


from  CK  U   compression  and  extension  tests.   No  correction  was 
made  to  the  field  vane  values  and  strength  anisotropy  was  con- 
sidered by  using  the  average  values  of  shear  strength  in  the 
horizontal  and  vertical  direction.   The  required  movements  were 
predicted  assuming  that  the  deformations  resulting  from  the  ad- 
ditional fill  would  occur  at  essentially  constant  volume.   Two 
different  methods  of  analyses  were  considered: 

(1)  Analysis  based  on  the  theory  of  elasticity 

(Boussinesq  solution)  and 

(2)  Analysis  using  the  finite  element  technique. 
Neither  method  is  considered  an  exact  solution.   In  the  elastic 
analysis,  a  value  (or  values)  must  be  assigned  to  the  modulus 
of  elasticity  (E)  and  Poisson's  ratio  (u) .   Poisson's  ratio  was 
taken  as  0.5  for  saturated  clay.   In  the  finite  element  approach, 
E,  u  and  yield  stress  were  selected  for  an  elastic  plastic  stress 
versus  strain  relationship.   Furnished  data  which  were  used  to 
arrive  at  pore  pressure  parameters  consisted  of  the  piezometer 


B  and  C  readings  and  the  CK  U   triaxial  compression  test  results. 
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Verne  C.  McGuffey  and  Raymond  L.  Gemme 
METHODS  EMPLOYED 

Additional  Height  of  Fill 

The  strength  of  the  clay  layer  used  in  the  stability 
analysis  was  developed  from  data  from  field  vane  tests  at 
Sta.  244  and  the  section  to  be  loaded  (Sta.  263)  and  from 
laboratory  CK  U  and  CIU  tests  run  on  samples  from  Sta.  263 
and  Sta.  24  6.   As  can  be  seen  in  Figure  IV  -  15/   vane  strengths 
at  Sta.  244  measured  in  1972  show  essentially  the  same  pro- 
file that  existed  in  1967  prior  to  construction.   It  was  assumed 
that  the  vane  strengths  at  Sta.  263  had  also  remained  essen- 
tially constant  with  time.   From  CIU  tests  run  with  a  consoli- 
dation pressure  four  times  the  maximum  past  pressure  (P  )  of 

P 

the  sample  a  shear  strength  ratio  C/P*  ,  for  the  normally  con- 
solidated clay  was  determined  to  be  0.28.   It  was  assumed  that 
this  ratio  was  the  same  for  CK  U  conditions  based  on  experience 
with  New  York  State  clays  where  a  maximum  difference  of  5%  has 
been  noted  between  strengths  measured  by  the  two  test  procedures . 
From  consolidation  test  data,  the  profile  of  maximum  past 

pressure  marked  P   (1967)  in  the  left  half  of  Figure  IV  -  16 
P 

was  developed  by  discounting  low  values  by  assuming  they  resulted 
from  sample  disturbance.   Since  the  clay  below  elevation  -60  ft 
is  thus  shown  to  be  slightly  overconsolidated,  a  shear  strength 
ratio  of  0.33  was  used  to  develop  the  strength  curve  marked  CK  U 
in  the  right  half  of  Figure  IV  -  16.    This  value  is  the  same 
factor  obtained  at  the  centerline  of  Sta.  246  in  the  slightly 
overconsolidated  range. 


-32- 


The  appropriate  strength  data  for  the  overconsolidated 
region  above  elevation  -60  feet  was  developed  using  the  pre- 
compression  ratio  (P  /P  )  with  Figure  IV-17  to  determine  the 

ir         *-' 

strength  ratio  (S  /C  )  of  the  normally  consolidated  strength 

to  the  overconsolidated  strength  at  the  same  effective  stress, 

This  strength  ratio  (S  /C  )  is  then  multiplied  by  the  shear 
u  u 

strength  ratio  (C/P  )  and  the  existing  effective  stress  to  de- 
termine the  strength.   Figure  IV-17  shows  a  band  of  values 
which  has  been  developed  with  experience  for  New  York  State 
clays.   The  precompression  ratio-strength  ratio  data  from  Sta, 
2  4  6  is  plotted  within  the  band  shown . 

It  was  concluded  that  the  vane  strengths  could  remain 
essentially  constant  due  to  a  cancellation  of  strength  gains 
due  to  consolidation  and  strength  losses  due  to  long  term  un- 
drained  strains  and  that  since  the  clay  is  in  a  state  of  plas- 
tic flow,  the  CK  U  strength  would  be  too  high.   Stability  an- 
alyses were  performed  using  three  different  assumptions  for 
the  strength  of  the  clay  layer: 

1.  Use  the  vane  strengths  (Figure  IV-16)  with  no 
allowance  for  strength  gain. 

2.  Use  the  average  of  the  vane  and  CK  U  strengths 
as  shown  in  Figure  IV-16  with  no  allowance  for 
strength  gain. 

3.  Use  the  average  of  the  vane  and  CK  U  strengths 
adjusted  to  reflect  the  strength  gain  due  to 
consolidation. 
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The  analyses  were  performed  using  a  Modified  Swedish  Wedge 
Analysis.   Figure  IV-18  shows  the  critical  wedge  for  Case  1 
noted  above.   A  circular  arch  analysis  was  also  run  using 
Bishop's  Method  which  yielded  approximately  the  same  results. 
Movements  Due  to  Adding  Six  Feet  of  Fill 
The  predictions  of  horizontal  movement  of  slope  indica- 
tors SI-3  and  SI-4,  the  settlement  of  SP-1  and  the  heave  of 
heave  points  H-l  and  H-2  were  all  derived  through  interpreta- 
tion of  the  data  from  the  filling  of  the  last  six  feet  of  fill 
placed  in  1969.   The  movements  of  the  slope  indicators  at  the 
required  points  and  the  settlement  plate  due  to  the  previous 
filling  are  shown  in  Table  IV-7.   Since  the  change  in  load  to 
the  foundation  due  to  the  additional  loading  would  be  less  than 
the  previous  loading  due  to  lower  unit  weight  of  the  new  fill 
material  and  a  narrower  crest,  the  measured  deformations  were 
divided  by  a  factor  of  1.18  (39/33)  to  account  for  the  lower 
load.   Similarly,  a  correction  factor  (0.9)  was  applied  to  account 
for  larger  deformations  due  to  the  increased  rate  of  loading 
anticipated.   The  horizontal  deformations  and  settlement  arrived 
at  from  applying  this  procedure  are  also  shown  in  Table  IV-7. 

Since  no  heave  measurements  were  made  during  the  previous 
loading,  predictions  of  heave  were  made  using  the  horizontal  de- 
formations of  slope  indicators  SI-3  and  SI-4  near  the  surface. 
Because  the  Bishop's  circle  gave  a  consistent  result  for  stabil- 
ity, it  was  used  to  extrapolate  the  horizontal  deformations  to 
estimate  the  heave  as  shown  in  Figure  IV-19 .   From  the  center  of 
the   critical  circle  near  surface  slope  indicator  points  were 
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found  so  that  the  rotation  at  SI-3  and  SI-4  would  be  the  same 
to  produce  the  measured  horizontal  movement.   The  vertical 
movement  (heave)  associated  with  that  rotation  was  then  plotted 
together  with  the  assumption  of  no  heave  below  the  critical 
center  and  no  movement  in  the  canal  to  derive  the  heave  vs.  dis- 
tance plot  shown  at  the  bottom  of  Figure  IV-19.   The  predicted 
heaves  for  points  H-l  and  H-2  were  then  determined  from  the 
derived  plot. 

Additional  Pore  Pressures 

The  pore  pressures  were  predicted  based  on  analysis  of 
the  pore  pressure  response  beneath  the  center-line  during  the 
original  loading.   Figure  IV-20  shows  the  piezometer  data  ex- 
trapolating back  to  the  anticipated  pore  pressure  under  instan- 
taneous loading.   It  was  felt  that  the  data  from  Piezometer  A 
were  probably  in  error  on  the  low  side  and  since  Piezometer  C 
was  damaged,  the  data  from  Piezometer  B  were  used  in  analyses. 
Table  IV-8  shows  the  change  in  Piezometer  B  which  would  have 
occurred  due  to  instantaneous  loading  for  each  of  the  three  pre- 
vious loading  periods  and  the  increase  in  load  with  two  assumed 
values  of  fill  density.   From  this  procedure,  it  was  then  con- 
cluded that  the  pore  pressure  response  as  measured  in  the  pie- 
zometers would  be  approximately  equal  to  the  increase  in  total 
stress  for  an  assumed  instantaneous  loading.   Since  the  loading 
would  be  much  faster  than  in  the  previous  loadings,  pore  pres- 
sure predictions  for  both  the  additional  6  feet  of  fill  and  at 
failure  were  made  on  the  assumption  that  the  increase  in  pore 
pressure  would  equal  the  increase  in  effective  stress. 
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PREDICTED  QUANTITIES 

Additional  Height  of  Fill 

From  the  stability  analysis  outlined  above/  a  plot  of 
additional  fill  height  versus  safety  factor  was  determined  for 
each  of  the  assumptions  with  regard  to  shear  strength  as  shown 
in  Figure  IV-21.   After  consideration  of  the  failure  mechanism 
and  the  vane  strength  data  at  Sta.  244,  it  was  decided  that  the 
field  vane  data  would  give  the  best  prediction.   Accordingly, 
the  height  of  fill  to  cause  failure  was  predicted  as  13.5  feet. 

Movements  Due  to  Adding  6  Feet  of  Fill 

The  calculated  movements  of  the  slope  indicator  points 
shown  in  Table  IV- 7  were  based  on  the  assumption  that  there 
would  be  no  structural  improvement  of  the  foundation  soils  asso- 
ciated with  a  strength  gain.   However,  in  estimating  the  defor- 
mations, it  was  felt  that  some  improvement  would  take  place  at 
least  for  Slope  Indicator,  SI-3.   Accordingly,  the  calculated 
values  were  modified  as  follows  to  give  the  predicted  deforma- 
tions : 


Cal 

culated 

Prediction 

SI-3 

(-30') 

1" 

0.25" 

SI-3 

(-70') 

1" 

0.25" 

SI-4 

(-30') 

0.25" 

0.25" 

SI-4 

(-70') 

0 

0 

Similarly,  the  calculated  settlement  of  3.5"  was  decreased  to  a 
predicted  settlement  of  3"  due  to  some  structural  improvement. 
The  predicted  heave  quantities  at  points  H-l  and  H-2  were  rounded 
to  1.0". 
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Additional  Pore  Pressures 

As  noted  above,  the  pore  pressures  were  assumed  to  be 
equal  to  the  change  in  total  stress: 


u  (ft) 

at 

6 

'  of 

fill 

at 

failure  (13.5') 

P-3 

7.3 

14.8 

P-4 

7.2 

13.9 

P-5 

4.9 

7.6 

DISCUSSION 

The  presentation  was  divided  into  two  parts  with  Mr. 
McGuffey  presenting  generalized  information  on  the  embankment 
design  process  and  the  practice  of  the  Bureau  of  Soil  Mechanics 
of  the  New  York  Department  of  Transportation  and  Mr.  Gemme  pre- 
senting the  prediction  methods  and  results.   In  Mr.  McGuffey 's 
presentation,  he  stressed  the  importance  of  experience  and  the 
prediction  was  influenced  substantially  by  the  predictors'  ex- 
perience with  Albany  Clay  which  was  reported  to  be  similar  to 
Boston  Blue  Clay.   In  discussing  the  design  process,  he  noted 
that  analysis  is  only  one  of  several  factors  contributing  to  the 
final  predictions  and  design  and  by  no  means  the  largest  portion 
of  the  effort  nor  the  largest  source  of  error.   Among  the  other 
factors  in  the  design  process,  he  mentioned  the  recent  history 
of  the  site  and  noted  the  possibility  that  the  foundation  was 
more  compressible  than  the  section  would  indicate   due  to  incom- 
plete  removal  of  the  peat  during  foundation  preparation.   How- 
ever, no  apparent  consideration  was  given  to  this  point  in  pre- 
paring the  centerline  settlement  prediction. 

The  consideration  of  strength  gain  due  to  consolidation  of 
the  embankment  was  a  critical  factor  in  the  prediction  made  (see 
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Figure  IV-21) .   It  was  not  considered  in  determining  the  height 
of  fill  and  consequently  the  pore  pressures  at  failure  nor  for 
the  estimate  of  heave.   It  was,  however,  considered  in  the  case 
of  lateral  movements  and  to  a  lesser  extent  in  the  settlement 
of  SP-1.   From  the  measured  results,  it  is  apparent  that  a  con- 
sistent application  of  a  gain  in  strength  and  structural  improve- 
ment due  to  consolidation  effects  would  have  improved  (strength 
gain)  or  weakened  (no  strength  gain)  the  prediction.   The  McGuffey- 
Gemme  predictions  are  summarized  in  Table  IV- 9 . 
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J.  M.  Duncan 

Prof.  Duncan  based  his  deformation  predictions  on  finite 
element  computations,  his  stability  prediction  on  limiting  equi- 
librium methods  and  his  pore  pressure  predictions  on  evaluation 
of  the  field  data.   Professor  Duncan's  predictions  are  summarized 
in  Table  IV-10. 

Consolidation  of  the  Clay 

Initially,  computations  were  made  to  determine  the  percent- 
age of  primary  consolidation  prior  to  new  construction  and  using 
the  results  of  these  calculations,  the  effects  of  undrained 
strength  and  deformation  properties  of  the  clay  were  estimated.   Set- 
tlement calculations  showed  that  at  most  (assuming  full  consoli- 
dation) only  1.6  ft  of  the  observed  3.8  ft  of  settlement  could 
be  accounted  for  by  consolidation.   The  remainder  was  undrained 
movement,  much  of  which  was  time -dependent.   Creep  effects,  there- 
fore, were  thought  to  be  very  important  for  Boston  Blue  Clay. 
Because  of  the  large  amount  of  settlement  due  to  undrained  creep, 
the  measured  pore  pressures  were  used  to  estimate  the  degree  of 
consolidation.   By  varying  c   (as  a  function  of  OCR) ,  the  ele- 
vation of  the  static  water  table ,  and  the  initial  excess  pore 
pressure,  good  agreement  between  computed  pore  pressures  and 
observed  pore  pressures  (for  July,  1974)  was  obtained.   The  con- 
clusions from  this  matching  were: 

1.  cv  =  35  (OCR)1,5  ft2/yr 

2.  Au        =  1.25   Aa, 

excess  1 
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3.  The  effective, average,  static,  water  table  elevation 
is  +2.5  ft  MSL  , 

4.  Of  the  3.8  ft  observed  settlement,  0.5  ft  was  due 
to  primary  consolidation. 

Creep 

Since  no  published  data  for  creep  tests  on  Boston  Blue  Clay 
were  available,  Professor  Duncan  requested  undisturbed  samples 
so  that  he  could  conduct  his  own  creep  tests.   Two  CIU   tests 
were  performed:   (a)  the  first  was  loaded  to  a  deviator  stress 
corresponding  to  40%  of  the  short-term  strength,  and  (b)  the 
second  was  loaded  to  a  deviator  stress  corresponding  to  60%  of 
the  short-term  strength.   These  loads  were  maintained  for  three 
weeks  after  which  the  observed  strains   were  2.5  to  3  times  the 
measured  strains  at  the  same  stress  levels  in  a  quick  test. 
To  account  for  these  larger  strains,  Professor  Duncan  reduced 
the  value  of  the  initial  modulus,  E.,  used  in  the  analyses.   For 
a  quick  CIU  test,  the  observed  value  of  E./s  was  approximately 
800.   The  one-day  creep  test  gave  E./s  -  400  while  three  weeks 
of  creep  resulted  in  E./s   -  250.   For  the  finite  element  anal- 

r  1/   U 

yses,  Professor  Duncan  used  E./s   =  400.   The  value  of  s   used 

1   u  u 

in  the  calculations  was  an  average  based  on  relationships  given 
in  Bulletin  3  (after  Davis  and  Christian,  1971) . 

Finite  Element  Model 

A  symmetrical  mesh,  Figure  IV-22,  was  constructed  for  pre- 
liminary work.   The  original  intention  was  to  use  a  complete 
non-symmetrical  mesh  later.   Close  agreement  between  observed 
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and  calculated  results  (Table  IV-8)  led  to  use  of  a  simple 
correction  procedure  and  the  symmetrical  mesh  instead  of  a 
larger,  unsymmetrical  mesh.   In  the  finite  element  analyses, 
the  hyperbolic  stress-strain  model  was  used  for  both  the  clay 
and  the  sand.   The  clay  was  characterized  by  values  of  E . /s  = 
400  and  s  /p  =  0 . 7  s   /p.   For  the  sand,  the  test  data  pro- 
vided by  MIT  were  inconsistent  and  estimated  values  based  on 
data  for  similar  sand  were  used  in  the  calculations.   The  same 
properties  were  used  for  the  natural  sand,  the  sand  working 
mat  and  the  embankment  sand.   Very  low  stiffness  was  assigned 
to  the  peat. 

The  predicted  values  of  settlement  and  horizontal  movement 
due  to  6  ft  of  additional  fill  were  calculated  using  values  of 
clay  strength  corrected  for  the  increased  effective  stress  re- 
sulting from  consolidation  under  the  weight  of  the  embankment. 
Deformations  were  computed  using  the  finite  element  method  and 
correcting  the  computed  values  with  factors  obtained  from  earlier 
comparisons.   Values  for  additional  pore  pressures  were  computed 
using  the  relationship  Au  =  1.25  Aa,  with  the  Aa,  taken  from 
the  finite  element  analyses. 

Stability 

Professor  Duncan  expressed  the  opinion  that  the  finite  ele- 
ment method,  with  hyperbolic,  elastic  stress-strain  relationships 
cannot  be  expected  to  give  reliable  results  for  conditions  near 
failure.   The  hyperbolic  stress-strain  relationship  is  useful 
for  calculating  movements  at  moderate  stress  levels  but  is  not  a 

good  method  of  predicting  stability  and  failure  loads.   The 
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stability  computations  for  his  prediction  of  failure  height 
were  made  using  the  Modified  Bishop  Method  of  analysis  with 
undrained  strengths  for  the  clay  and  drained  strengths  for 
the  sand. 

To  determine  how  creep  affected  the  ultimate  strength  of 
the  clay,  both  the  40%  and  60%  CIU  tests  were  loaded  to  90%  of 
short-term  strength  and  both  failed  within  1  1/2  days.   The 
strengths  used  in  the  stability  calculations  were  the  same 
averaged  s   values  used  in  the  finite  element  calculations  but 
reduced  by  10%  to  account  for  loss  of  strength  due  to  creep. 
These  reduced  strengths  were  not  used  in  the  determination  of 
E.  for  the  finite  element  analyses  since  the  creep  test  strengths 
were  not  available  when  the  finite  element  analyses  were  performed 
The  calculated  values  of  factor  of  safety  were  plotted  against 
fill  height  for  both  full  s   and  90%  s   accounting  for  10%  creep 
strength  loss  (Figure  IV-2  3) .   The  predicted  value  of  embankment 
height  at  failure  was  chosen  from  the  90%  s   curve  at  a  factor 
of  safety  equal  to  1.0,  and  then  discounted  from  14  ft  to  12  ft 
of  additional  fill  because  Professor  Duncan  suspected  that  the 
laboratory  tests  might  result  in  strengths  which  were  higher  than 
the  field  strengths. 

Failure  Pore  Pressures 

After  prefacing  his  remarks  with  a  disclaimer  on  the  use- 
fulness of  pore  pesssures  in  an  undrained  analysis,  Professor 
Duncan  presented  his  failure  pore  pressures.   They  were  computed 
using  the  Au  =  1.25  Aa   relationship  obtained  from  the  consolida- 
tion computations.   The  change  in  maximum  principal  stress  was 
computed  using  an  additional  fill  height  of  12  ft. 
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Summary 

Professor  Duncan's  major  concerns  were  how  to  account 
properly  for  creep  effects  in  Boston  Blue  Clay.   Using  the 
results  of  the  two  creep  tests  he  had  performed  as  a  guide, 
he  reduced  the  slope  of  the  stress-strain  curve  and  the  value 
of  s  .   His  predictions  were  based  on  the  finite  element 
method  (modified  on  the  basis  of  comparisons  between  measured 
and  calculated  movements) ,  an  evaluation  of  field  data  (checked 
with  consolidation  computations  for  pore  pressures)  and  limit- 
ing equilibrium  analyses  of  stability. 

Appendix  D  presents  a  more  detailed  description  of  Pro- 
fessor Duncan's  predictions. 
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R.  C.  Kirby 

The  approach  used  by  this  team,  Woodward-Moorhouse  and 
Associates,  Inc.  (WMAI)  was  that  of  a  consulting  engineer,  con- 
sidering the  researchers  to  be  the  client.   The  clients'  "needs" 
were  predictions  of  deformations,  pore  pressures  and  the  failure 
load  on  the  clay.   A  detailed  treatment  of  WMAI * s  prediction 
is  found  in  Appendix  E. 

It  is  likely  that  a  consultant  would  make  the  prediction 
of  deformation  and  pore  pressures  due  to  the  addition  of  six 
feet  of  fill  by  first  replotting  the  available  field  measure- 
ments data  and  then  extrapolating  them  (using  judgment) .   The 
prediction  of  the  height  of  fill  to  cause  failure  would  be  con- 
sidered as  a  practical  tool  only  if  extrapolating  the  measured 
data  was  considered  very  uncertain,  and  the  consequences  of  be- 
ing wrong  included  a  severe  cost  penalty  or  an  unacceptable  risk 
to  the  client. 

In  order  to  arrive  at  the  needed  predictions  and  also  evalu- 
ate methods  of  prediction  performances,  the  project  approach 
shown  on  Fig.  IV-24  was  adopted.   The  key  characteristics  of 
this  approach  is  an  interim  evaluation  of  the  Finite  Element 
Method  by  determining  if  it  could  successfully  hindcast  the  mea- 
sured performance  of  the  embankment  previously  constructed  to  elev. 
40.   This  project  approach  was  developed  by  the  project  engineer 
(Mr.  Kirby)  and  reviewed  and  approved  by  senior  level  members 
of  WMAI.   The  judgments  required  at  the  decision  points  were 
made  by  the  project  engineer;  these  judgments  were  also  reviewed. 
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If  the  Finite  Element  Method  proved  satisfactory  for  hind- 
casting  performance,  it  would  be  used  as  a  prediction  tool  in 
combination  with  other  methods.   If  the  Finite  Element  Method 
proved  unsatisfactory,  predictions  would  be  made  by  extrapo- 
lating the  measured  data  and  the  use  of  limiting  equilibrium 
methods.   Although  not  shown  on  the  project  approach  diagram, 
"quick  and  dirty"  checks  of  the  analyses  were  included  whenever 
possible. 

WMAI  was  unable  to  hindcast  the  measured  performance  of 
the  embankment  constructed  to  elev.  40  using  the  finite  element 
method.   The  details  of  the  finite  element  analysis  are  included 
in  Appendix  E.   The  predictions  are  summarized  in  Table  IV-11. 

Movements  and  pore  pressures  were  predicted  by  first  inter- 
preting and  then  extrapolating  the  measured  settlement,  horizontal 
movement  and  pore  pressure. 

Construction  from  elev.  10  to  elev.  34  required  approximately 
five  months  to  complete  and  a  portion  of  the  settlement  measured 
during  construction  is  due  to  consolidation  of  the  clay.   On  the 
other  hand,  the  addition  of  six  feet  of  fill  in  September  19  74 
proceeded  at  a  rate  approximately  five  times  faster  than  the  con- 
struction from  elev.  +10  to  elev.  +34.   Consequently,  there  was 
a  need  to  separate  the  undrained  or  immediate  settlement  from  the 
total  settlement.   A  rate  of  settlement  analysis  was  made  using 
the  computer  program  PROGRS  (Schiffman  and  Stein,  1969)  in  order 
to  determine  what  the  shape  of  a  theoretical  settlement  vs.  square 
root  of  elapsed  time  plot  would  be.   The  rate  analysis  included 
the  effects  of  loading  history,  variation  of  cv  with  depth,  and 
approximate  corrections  for  two-dimensional  strain  and  drainage. 
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Both  piezometers  (where  predictions  were  required)  were  approx- 
imately 20  feet  from  the  nearest  drainage  layer  and  dissipation 
of  pore  pressure  during  construction  at  these  locations  was  con- 
sidered negligible.   The  measured  horizontal  movements  during 
construction  were  probably  affected  somewhat  by  consolidation 
of  the  clay  during  construction.   However,  no  reasonably  simple 
means  of  investigating  this  effect  are  available  and  no  correc- 
tions were  made. 

Recognizing  the  uncertainty  of  any  prediction,  it  was  de- 
cided to  estimate  the  deformation  caused  by  adding  six  feet  of 
fill  in  1974  to  be  the  same  as  the  measured  deformations  due  to 
raising  the  embankment  from  elev.  34  to  elev.  40  in  1968.   The 
predicted  additional  pore  pressure  at  piezometers  P-3  and  P-4 
are  based  on  extrapolating  the  measured  pore  pressures  at  pie- 
zometers A  and  B.   These  data   show     a  two-feet  increase  in 
total  head  for  each  one-foot  increase  in  fill  elevation.   For  the 
additional  6  ft  of  fill,  a  12-foot  increase  in  total  head  at 
both  piezometers  P-3  and  P-4  was  predicted.   The  additional  pore 
pressure  at  P-6  (which  is  located  near  the  toe  of  the  embankment) 
could  not  be  predicted  by  extrapolating  measured  data  during  con- 
struction since  this  data  does  not  exist.   The  additional  pore 
pressure  was  estimated  analytically  from  the  change  in  average 
total  normal  stress  at  that  location  which  was  determined  from 
the  finite  element  analysis.   The  same  logic  described  above  for 
extrapolating  the  settlement  data  was  applied  to  extrapolating 
the  horizontal  movement  data.   No  prediction  was  made  for  the 
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heave  at  H-l  and  H-2.   This  prediction  can  only  be  made  analy- 
tically because  no  previous  measurements  are  available.   WMAI * s 
analytical  model  (the  finite  element  method)  proved  unreliable 
for  predicting  deformation,  and  other  simpler  analytical  methods 
were  deemed  inappropriate. 

The  plane  strain  active  and  passive  strengths  are  con- 
sidered a  best  estimate  for  the  preconstruction  undrained  strength 
but  have  not  been  verified  by  full-scale  local  tests  to  failure. 
Analyses  were  made  using  the  ICES  program  LEASE  (Bailey  and 
Christian,  1969).   The  program  considers   circular  failure  arcs 
and  computes  the  factor  of  safety  using  both  the  Fellenius  and 
Simplified  Bishop  method  of  slices. 

Effective  stress  parameters  were  input  for  the  embankment 
fill  and  natural  sand  (c  =  0,  <j>  =  40°)  and  undrained  shear  strengths 
were  used  for  the  clay.   The  analyses  were  made  for  three  dif- 
ferent profiles  of  undrained  shear  strength:  Case  1- -field  vane 
strengths;  Case  2--preconstruction  plane  strain  active  and  passive 
strengths,  and  Case  3 --plane  strain  active  and  passive  strength 
increased  to  approximately  account  for  the  increase  in  undrained 
shear  strength  due  to  consolidation. 

By  analyzing  the  performance  of  the  embankment  constructed 
in  196  8  and  1969  to  elev.  40  using  the  finite  element  method, 
WMAI  could  not  hindcast  the  measured  behavior  without  making 
major  revisions  in  the  input  parameters.   The  input  parameters 
chosen  for  the  finite  element  analyses  are  considred  best  esti- 
mates and  major  changes  in  these  parameters  are  not  justified. 
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Moreover,  because  of  the  many  parameters  which  could  be  changed 
to  give  good  agreement  between  measured  and  predicted  perfor- 
mance, the  significance  of  such  agreement  would  be  questionable, 
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G .  P .  Raymond 

The  prediction  of  performance  of  the  embankment  founded 
on  Boston  Blue  Clay  made  by  Dr.  Raymond  was  essentially  based 
on  a  numerical  model:   a  finite  element  program.   The  strength 
properties  used  for  the  undrained  analysis  were  obtained  from 
field  vane  shear  strengths  with  no  reduction  to  account  for  Pi . 
The  results  of  the  vane  tests,  together  with  the  results  of  un- 
confined  compression  and  UU  tests  are  presented  in  Figures  IV- 
25  and  IV-26.   Fig.  IV-26,  in  addition,  contrasts  the  field  vane 
strength  prior  to  any  construction  to  the  strengths  after  Stage 
1  was  built.   This  comparison  seems  to  indicate  no  increase  in 
shear  strength  with  time  after  the  construction  of  Stage  1. 
The  straight  line  on  both  figures  represents  the  assumed  soil 
strength.   Past  experience  with  Canadian  clays  has  suggested  that 
the  results  of  laboratory  tests  performed  on  undisturbed  samples 
consolidated  to  pressures  higher  than  the  maximum  preconsolida- 
tion  pressures  are  not  indicative  of  the  soil  in  situ  properties. 
Consolidating  to  higher  pressures  than  a    breaks  the  soil  struc- 
tures and  changes  the  soil  properties.   Strength  parameters  from 
field  vane  tests  on  Canadian  clays  are  more  representative  of 
real  in  situ  properties  than  results  based  on  laboratory  tests 
that  follow  the  technique  described  above. 

The  stress-strain  properties  needed  for  the  lateral  move- 
ment predictions  performed  by  the  finite  element  program  were  ob- 

al"a3 
tained  from  the  New  Liskeard  ( )  vs.  strain  data  and  the  MIT 
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vane  results.   Because  of  time  considerations,  some  of  the 
laboratory  tests  requested  could  not  be  performed.   The  tests 
requested  were  undrained  test,  consolidated  isotropically  to 
1000  psf,  sheared  to  4000  psf,  expanded  to  3/4  of  the  effective 
overburden  and  finally  sheared  well  beyond  failure  (to  about 
20%  strain) .   Data  from  these  tests  was  also  needed  for  pore 
pressure  predictions. 

The  New  Liskeard  failure  index  vs.  strain  is  presented 
in  Fig.  IV-2  7.   Boston  Blue  Clay  was  suspected  to  have  higher 
strength  than  the  New  Liskeard  varved  clay.  For  this  reason,  the  New 
Liskeard  strains  were  divided  by  a  suitable  factor  to  fit  the  MIT  data. 
This  was  accomplished  with  a  factor  of  E,    >  =  E  ,   ,  /(1.5)  . 

The  strength  parameters  of  the  fill  used  were  0  =  37   and 
Y   =12  0  pcf.   A  similar  analysis  on  an  embarkment  built  near 
Kars,  Ontario,  suggested  that  the  embarkment  stiffness  should 
be  neglected  in  the  stability  analysis.   This  approach  was  fol- 
lowed . 

According  to  the  finite  element  analysis  mode  for  6  ft  of 
additional  fill,  all  the  piezometers  had  either  failed  or  were 
very  close  to  failure.   To  predict  the  increase  in  pore  pressure 
a  value  of  A  equal  to  1.0  together  with  Skempton's  pore  pressure 
equations  were  used  for  all  the  piezometers.  In  addition  to 
this,  for  those   piezometers   that  had  already  yielded,  a  range 
of  predicted  pore  pressure  was  obtained  by  computing  an  upper 
value  based  on  an  empirical  equation.   The  equation  used  is  as 
follows : 
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A  U   =  a   vo   +  Aav   ~   3Su   -   U0 

Y  w 

where   a  =  effective  vertical  stress  prior  to  the 
addition  of  6 '  of  fill 

A  av   =  increment  in  total  vertical  load 

Su   =  undrained  strength  as  computed  using  field 
vane. 

w  =  unit  weight  of  water 

Uo  =  total  pore  pressure  prior  to  the  increment 
in  load . 

The  same  approach  was  used  for  predicting  pore  pressure  at 
failure. 

Conclusions: 

I.   Additional  fill  to  cause  failure: 

a.  Model  used 

1.   numerical  model  based  on  finite  element  approach 

b.  Strength  Parameter 

1.   Obtained  from  field  vane  test  with  no  correction 
to  account  for  PI 

c.  Embankment  stiffness  neglected 

d.  Prediction:   8  feet  of  additional  fill. 
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II.   Additional  Movement  for  6  feet  of  fill 

A.  Model  used 

1.   finite  element  program 

B.  stress-strain  data 
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1.   Obtained  from  New  Liskeard   , r     vs.  E 

together  with  the  vane  results. 


2 .   Correction  on  NL  strains  was  made  to  account 
BBC  higher  strength. 

C.  Prediction  of  lateral  movements: 

1.  SI  -  3  at   -  30  0.72  ft. 

2.  SI  -  3  at   -  70  0.72  ft. 

3.  SI  -  4  at   -  30  0.44  ft. 

4.  SI  -  4  at   -  70  0.00  ft. 

D.  Predictions  of  movements  of  surface 

1.  At  center  line        1.14  ft. 

2.  At  heave  pt.  H-l       0.40  ft. 

3.  At  heave  pt.  h-2       0.52  ft. 

III.   Additional  Pore  Pressure  for  6  ft  of  fill 

A.  Model 

1.   Skempton  pore  pressure  equation — lower  value 

3su 
2-   AV=  °vo  +  AV  T^T  -  Uo        —upper  value 

B.  Parameter 

1.  A  =  1.0  (assumed) 

2.  B  =  1.0  (saturated  soil) 

3.  su   =  field  vane  test 

4.  Stresses  and  existing  pore  pressure  as  predicted 
by  finite  element  program. 
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C.      Prediction   for    6    ft      of    fill1 


1.  P-3 

2.  P-4 

3.  P-6 


between 
between 
between 


Prediction  at  failure: 

1.  P-3         between 

2.  P-4         between 

3.  P-6         between 


10 '  -  43*  of  water 

9'  -  35'  of  water 

11'  -  of  water 

12'  -  46'  of  water 

11'  -  38'  of  water 

13'  of  water 


Additional  figures  and  data  used  in  Dr.  Raymond's  predic- 
tion  are   presented  in  Appendix  F. 

Professor  Raymond's  predictions  are  summarized  in  Table 
IV  -  12. 
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C.  P.  Wroth 

Dr.  Wroth  based  his  predictions  on  the  Cam -clay  model  of 
clay  behavior  developed  at  Cambridge  University. 

Information  on  the  Cam  -clay  model  of  soil  behavior  and 
details  of  the  prediction  are  presented  in  Appendix  G.   Para- 
meters for  the  Cam-clay  model  were  selected  from  laboratory 
tests  given  in  the  symposium  bulletins.   Values  used  in  the 
computations  to  model  Boston  Blue  Clay  were  Cc  =  0.3  38,  Cr  = 
0.138,  r=  2.15,  <j>'=  26.5°  and  v=  0.35  (see  Appendix  G  for  no- 
tation. ) 

Behavior  of  the  sand  fill  was  based  on  a  model  developed 
by  Simpson  (1973)  which  incorporates  provisions  for  strain  soften- 
ing and  dilation.   The  existing  model  with  parameters  chosen  to 
match  Leighton  Buzzard  sand  was  used  without  modification. 

The  finite  element  program  used  for  the  predictions  has 
many  special  features,  the  most  notable  being: 

1.  The  plane  strain  elements  used  are  simple  3-node 
constant  strain  triangles. 

2.  The  program  can  handle  strain-softening  and  dilation. 

3.  The  program  treats  soil  as  a  two-phase  material;  the 
behavior  of  the  soil  is  always  controlled  by  the 
effective  stresses  and  equilibrium  by  total  stres- 
ses.  The  pore-water  is  ascribed  an  appropriately 
large  (but  finite)  bulk  modulus,  so  that  for  the 
so-called  undrained  condition  there  has  to  be  compa- 
tibility between  the  small  volume  change  of  the  water 
and  that  of  the  soil  skeleton.   The  program  keeps 
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track  of  total  and  effective  stresses  for  every 
element.   For  example,  Fig.  IV-28  shows  the  effec- 
tive and  total  stress  paths  for  the  indicated  ele- 
ment positioned  under  the  center  line  of  the  em- 
bankment; the  stage  at  which  this  element  yielded 
is  clearly  shown  by  the  abrupt  change  in  the  direc- 
tion of  the  effective  stress  path  at  Y.   At  Z,  the 
element  has  reached  failure  and  is  being  held  in 
equilibrium  by  the  surrounding  elements  but  unable 
to  sustain  any  greater  shear  stress.   The  first 
part  XY  of  the  effective  stress  path  is  the  partly 
drained  simulation  of  the  combined  effects  of  the 
initial  construction  and  subsequent  five  years  of 
partial  consolidation;  the  method  is  described  below. 

4.  The  program  allows  for  the  effect  of  independent 
rotation  of  the  principal  axes  of  stress,  stress- 
increment  and  strain-increment. 

5.  The  program  includes  a  facility  for  refining  the 
mesh  of  elements  during  a  computation  so  that  pro- 
gressive failure  and  growth  of  failure  zones  may 
be  modelled. 

Matching  Field  Consolidation 

The  predictions  are  based  on  two  separate  computations.   The 
first  computation  was  modelling  the  construction  of  the  initial 

embankment  up  to  the  level  of  +40  ft  as  a  partially  drained  situa- 

5        2 
tion,  for  which  a  bulk  modulus  for  the  water  K   of  1  x  10   lbf/ft 
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was  adopted;  this  gave  a  reasonable  match  to  the  field  data 
available.   The  criterion  for  acceptability  was  when  adequate 
agreement  had  been  obtained  for  the  computed  horizontal  move- 
ments with  the  field  data  from  the  inclinometers  after  600 
days  (Aug.  1969).   The  computed  profiles  are  given  in  Fig.  IV-29. 

The  field  data  for  the  changes  that  occurred  between  600 
days  (Aug.  1969)  and  5  years  (Aug.  19  74)  showed  small  changes 
in  excess  pore  pressures  but  substantial  settlements  which  were 
attributed  mainly  to  undrained  creep.   Because  the  in  situ  ef- 
fective stresses  had  not  substantially  altered  during  this  per- 
iod, it  was  believed  that  the  computation  represented  as  well 
as  it  could  the  existing  stresses  in  the  ground  at  the  start  of 
the  test.   The  final  stage  was  to  model  the  filling  of  the  embank- 
ment up  to  the  predicted  failure  height  of  elevation  +61.   This 
was  carried  out  as  a  completely  undrained  condition  with  the  bulk 
modulus  of  water  of  2  x  106  lbf/ft  . 


Failure 

Failure  of  the  embankment  could  not  be  properly  modelled 
by  the  finite  element  program  because  the  Cam-clay  model  reaches 
failure  (as  defined  by  maximum  shear  stress)  for  normally  consol- 
idated soil  only  after  infinite  strain. 

Recourse  was  made  to  a  stability  analysis  using  the  classi- 
cal method  of  a  slip  circle.   A  limited  search  for  the  worst  circle 
was  made,  and  the  parameters  chosen  for  the  shear  strength  were 
(1)  fill  material  c'=0,  <t>'=30°,  Y=120  lbf/ft3;  (2)  sand  layer  in 
the  existing  ground  c'-0,  <j>'=35°,  Y=120  lbf/ft3;  (3)  Boston  Blue 
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2  2 

Clay  c   =8.8  lbf/in   (60  kN/m  ).   The  last  figure  was  chosen 

partly  on  the  quoted, in  situ  vane  shear  strengths,  and  partly 
on  the  results  of  undrained  pressuremeter  tests  carried  out 
at  the  site  by  Hughes  using  the  Camkometer  (self -boring  pres- 
suremeter) developed  at  Cambridge,  Wroth  and  Hughes  (1973) . 

No  allowance  was  made  for  anisotropy  of  strength  in  the 
clay;  a  10-ft  vertical  crack  was  assumed  at  the  top  of  the  fill 
where  the  slip  circle  reached  the  surface;  reduced  strength  was 
taken  for  the  small  length  of  arc  at  the  toe  of  the  slip  circle 
where  it  passed  through  the  peat. 

Dr.  Wroth 's  predictions  are  summarized  in  Table  IV-13 . 
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R.  H.  Bassett 

Dr.  Bassett' s  prediction  was  based  on  a  centrifuge  model 
of  the  1-95  embankment.   Appendix  H  presents  a  description  of 
the  centrifugal  model  technique. 

The  1-9  5  embankment  model  foundation  material  consisted 
of  powdered  Boston  Blue  Clay  supplied  by  MIT  and  remixed  with 
distilled  water  under  vacuum  at  1.7  to  1.9  times  the  liquid 
limit.   The  normal  consolidation  properties  of  this  remolded 
material  compared  well  with  the  reported  properties.   Study 
of  Fig.  15  in  the  Symposium  Bulletin  suggested  that  a  reason- 
able approximation  to  the  prototype  would  be  a  two-layer  model 

2 

with  OC  clay  with  a  Pc  of  5.8  kip/ft   from  -65  feet  downward. 

The  stiff  clay  from  -12'  to  -18'  elevation  (Pc  =  8.5  kip/ft2) 
was  omitted.  The  1-95  model  had  a  simplified  horizontal  till 
strata. 

Due  to  the  limited  capacity  of  the  available  centrifuge 
at  Luton,  England,  a  1/180  scale  model  was  tested  at  60  g's 
(i.e.,  the  model  represented  a  prototype  1/3  the  actual  size). 
It  was,  therefore,  necessary  to  consolidate  materials  at  stres- 
ses one-third  of  those  in  the  field.   The  slurry  was  K0  consoli- 
dated in  two  large  containers  (to  a  consolidation  stress  B^j  = 
0.62  ton/ft2  for  the  NC  clay  (Layer  B)  and  Xm  =  0.84  ton/ft2 
for  the  OC  clay  (Layer  A) )  and  then  suitable  layers  were  trimmed 
to  the  required  thickness  and  placed  in  the  test  package.   It 
was  later  found  that  the  actual  consolidation  pressures  were 
Bm  =  0.7  3  ton/ft2  and  A   =  1.15  ton/ft2.   The  error  in  the  lower 
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layer  was  very  little  but  in  the  upper  layer  it  was  consid- 
erable. 

The  model  was  run  in  the  centrifuge  at  60  g's  for  8  hours 
to  bring  it  to  stress  similarity  with  a  1/3  full  scale  pro- 
totype.  The  machine  was  stopped  and  an  embankment  construc- 
ted to  represent  the  original  work.   The  model  was  again  accel- 
erated to  represent  construction  as  shown  in  Figures  IV-30 
and  maintained  for  just  over  2  1/2  hours  to  represent  the  var- 
ious surcharge  loading  stages.   During  the  whole  sequence  of 
operations,  photographs  for  strain  measurements  and  pore  pres- 
sure transducer  readings  were  taken. 

Strain  measurements  for  the  construction  of  the  original 
embankment  show  some  indication  of  non-circular  characteristics. 
They  also  give  the  first  indication  of  multiple  failure  mecha- 
nisms acting  together.   At  the  end  of  the  original  embankment 
construction,  it  was  known  that  the  original  laboratory  consol- 
idation had  resulted  in  stiffer  materials  being  used  in  the  model 
than  the  ideal  although  the  relative  stiffnesses  would  be  of 
the  same  ratio.   In  order  to  obtain  a  scale  factor,  the  model 
and  prototype  performance  were  compared.   The  scale  comparison 
between  prototype  and  model  for  a  reasonable  fit  on  the  hori- 
zontal displacement  data  would  be  500  to  1  and  for  the  verti- 
cal displacement  1000  to  1,  not  180:1. 

The  embankment  surcharge  was  contained  in  a  hopper  opened 
by  a  hydraulic  piston  which  forms  part  of  the  test  box.   The 
hopper  was  actuated  for  three  seconds  about  130  minutes  after 
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the  original  construction  had  been  completed.   The  surcharge 
weight  turned  out  to  be  equivalent  to  8  feet  of  fill  at  25% 
moisture  content.   This  was  reasonably  close  to  the  requested 
6  foot  value.   The  displacement  results  in  terms  of  incremen- 
tal values  immediately  after  loading,  are  shown  in  Figures  IV- 
31  and  IV-32.   They  are  shown  as  actual  model  displacements 
with  no  allowances  for  stiffness  errors.   The  volumetric  and 
shear  strain  are  shown  in  Figures  IV-33  and  IV-34. 

Finally,  after  a  three-minute  pause  (equivalent  to  four 
weeks)  to  record  pore  water  pressures,  the  hopper  was  opened 
for  10  seconds  and  the  measurements  repeated.   The  load  was 
greater  than  intended  and  caused  considerable  widening  of  the 
base.   The  data  for  this  last  step  is  presented  in  Appendix  H. 

Dr.  Bassett's  predictions  are  summarized  in  Table  IV-14. 
The  displacement  predictions  were  obtained  by  multiplying  the 
actual  model  displacements  by  the  scaling  factor  obtained  be- 
fore from  comparison  with  actual  field  data. 
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TABLE  IV-2.   SUMMARY  OF  PREDICTIONS 
Poulos  &  Davis 


QUANTITY 

PREDICTED  VALUE 

Additional  height  of  fill  to 
cause  stability  failure 

27.0     ft 

Additional  horizontal  move- 
ments for  6  ft.  extra  fill 

SI-3    El.  -30 

El.  -70 
SI-4    El.  -30 

El.  -70 

0.17" 

0.016" 

0.47" 

o 

>   All 
eastwards 

Additional  settlement  of  SP-1 
for  6  ft  extra  fill 

0.76" 

Additional  pore  pressures  (in 
ft.  of  water)  for  6  ft.  extra 
fill 

P3 
P4 
P6 

8.8 
9.5 
8.1 

Additional  heave  for  6  ft 
extra  fill 

HI 
H2 

0.23" 
0.32" 

Additional  pore  pressures  at 
failure  (in  ft  of  water) 

P3 
P4 
P6 

17.9 
22.1 

12.1 
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TABLE   IV- 3.      W.W.    MOORE   -   PREDICTIONS 


QUANTITY 

PREDICTED  VALUE 

Additional  height  of  fill  to 
cause  stability  failure 

13    ft 

Additional  horizontal  move- 
ments for  6  ft.  extra  fill 

SI-3    El.  -30 

El.  -70 
SI-4    El.  -30 

El.  -70 

0.75")   All 

0. 75")  eastwards 

0.75") 

0.50") 

Additional  settlement  of  SP-1 
for  6  ft  extra  fill 

1.25" 

Additional  pore  pressures  (in 
ft.  of  water)  for  6  ft.  extra 
fill 

P3 
P4 
P6 

8.8 
7.2 
4.8 

Additional  heave  for  6  ft 
extra  fill 

HI 
H2 

1.0" 
0.5" 

Additional  pore  pressures  at 
failure  (in  ft  of  water) 

P3 
P4 
P6 

16.6 
16.6 
11.9 
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TABLE  IV-5.   C.C.  LADD  -  PREDICTIONS 


QUANTITY 

PREDICTED  VALUE 

Additional  height  of  fill  to 
cause  stability  failure 

23.0    ft. 

Additional  horizontal  move- 
ments for  6  ft.  extra  fill 

SI-3   El.  -30 

El.  -70 
SI-4   El.  -30 

El.  -70 

0.7"  ] 

0.35"     All 
0.45"  I  eastwards 
0     J 

Additional  settlement  of  SP-1 
for  6  ft.  extra  fill 

0.8" 

Additional  pore  pressures  (in 
ft.  of  water)  for  6  ft.  extra 
fill 

P3 
P4 
P6 

14 
14 
10 

Additional  heave  for  6  ft. 
extra  fill 

HI 
H2 

0.5"  ' 
0.4" 

Additional  pore  pressures  at 
failure  (in  ft.  of  water) 

P3 
P4 
P6 

37 
39 
25 
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TABLE   IV- 6.      Z.C.    MOH   -   PREDICTIONS 


QUANTITY 

PREDICTED  VALUE 

Additional  height  of  fill  to 
cause  stability  failure 

9 

ft. 

Additional  horizontal  move- 
ments for  6  ft.  extra  fill 

SI-3   El.  -30 

El.  -70 
SI-4   El.  -30 

El.  -70 

1.5"   1 

1.25" 

1.0" 

0   "   J 

All 
eastwards 

Additional  settlement  of  SP-1 
for  6  ft.  extra  fill 

3.0" 

Additional  pore  pressures  (in 
ft.  of  water)  for  6  ft.  extra 
fill 

P3 
P4 
P6 

9.5 
9.2 

11.5 

Additional  heave  for  6  ft. 
extra  fill 

HI 
H2 

0.25" 
1.25" 

Additional  pore  pressures  at 
failure  (in  ft.  of  water) 

P3 
P4 
P6 

12.8 
12.5 
17.1 
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TABLE  IV- 9.   McGUFFEY  &  GEMME  -  PREDICTIONS 


QUANTITY 

PREDICTED  VALUE 

Additional  height  of  fill  to 
cause  stability  failure 

13.5  ft 

Additional  horizontal  move- 
ments for  6  ft.  extra  fill 

SI-3   El.  -30 

El.  -70 
SI-4   El.  -30 

El.  -70 

.25"  1 

.25"     All 
. 25"  [  eastwards 
0     J 

Additional  settlement  of  SP-1 
for  6  ft.  extra  fill 

3.0" 

Additional  pore  pressures  (in 
ft.  of  water)  for  6  ft.  extra 
fill 

P3 
P4 
P6 

7.3 
7.2 
4.9 

Additional  heave  for  6  ft. 
extra  fill 

HI 
H2 

1.0" 
1.0" 

Additional  pore  pressures  at 
failure  (in  ft.  of  water) 

P3 
P4 
P6 

14.8 

13.9 

7.6 
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TABLE   IV- 10.      J.M.    DUNCAN   -   PREDICTIONS 


QUANTITY 

PREDICTED  VALUE 

Additional  height  of  fill  to 
cause  stability  failure 

12    ft. 

Additional  horizontal  move- 
ments for  6  ft.  extra  fill 

SI-3   El.  -30 

El.  -70 
SI-4   El.  -30 

El.  -70 

1.2"  | 

0.6"     All 
0.6"  1  eastwards 
0    j 

Additional  settlement  of  SP-1 
for  6  ft.  extra  fill 

1.8" 

Additional  pore  pressures  (in 
ft.  of  water)  for  6  ft.  extra 
fill 

P3 
P4 
P6 

19.22 
20.82 
12.81 

Additional  heave  for  6  ft, 
extra  fill 

HI 
H2 

0.6" 
1.2" 

Additional  pore  pressures  at 
failure  (in  ft.  of  water) 

P3 
P4 
P6 

38.44 
41.65 
25.63 
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TABLE  IV-11.   R.C.  KIRBY  -  PREDICTIONS 


QUANTITY 

PREDICTED  VALUE 

Additional  height  of  fill  to 
cause  stability  failure 

23.5  ft. 

Additional  horizontal  move- 
ments for  6  ft.  extra  fill 

SI-3   El.  -30 

El.  -70 
SI-4   El.  -30 

El.  -70 

-1  9"  ) 

L- Z       All 
-0  8"     Aii 
~X' 2n  r  eastwards 

Additional  settlement  of  SP-1 
for  6  ft.  extra  fill 

-5" 

Additional  pore  pressures  (in 
ft.  of  water)  for  6  ft.  extra 
fill 

P3 
P4 
P6 

12   ft. 

12   ft. 

2.5  ft. 

Additional  heave  for  6  ft. 
extra  fill 

HI 
H2 

No  Prediction 

Additional  pore  pressures  at 
failure  (in  ft.  of  water.) 

P3 
P4 
P6 

25.5  ft. 
25.5  ft. 
No  Prediction 
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TABLE  IV- 12.   G.P.  RAYMOND  -  PREDICTIONS 


QUANTITY 

PREDICTED  VALUE 

Additional  height  of  fill  to 
cause  stability  failure 

8   ft. 

Additional  horizontal  move- 
ments for  6  ft.  extra  fill 

SI-3   El.  -30 

El.  -70 
SI-4   El.  -30 

El.  -70 

8.6") 

8.6"     All 
5.3"  |  eastwards 
0    J 

Additional  settlement  of  SP-1 
for  6  ft.  extra  fill 

13.7" 

Additional  pore  pressures  (in 
ft.  of  water)  for  6  ft.  extra 
fill 

P3 
P4 
P6 

26 
22 
11 

Additional  heave  for  6  ft. 
extra  fill 

HI 

H2 

4.8" 
6.2" 

Additional  pore  pressures  at 
failure  (in  ft.  of  water) 

P3 
P4 
P6 

29 
24 
13 
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TABLE  IV-13.   C.P.  WROTH  -  PREDICTIONS 


QUANTITY 

PREDICTED  VALUE 

Additional  height  of  fill  to 
cause  stability  failure 

21   ft. 

Additional  horizontal  move- 
ments for  6  ft.  extra  fill 

SI-3   El.  -30 

El.  -70 
SI-4   El.  -30 

El.  -70 

1.6"  1 

0.6"     All 
0.7"  [  eastwards 
0    J 

Additional  settlement  of  SP-1 
for  6  ft.  extra  fill 

2.6" 

Additional  pore  pressures  (in 
ft.  of  water)  for  6  ft.  extra 
fill 

P3 
P4 
P6 

9 

9.6 

5.8 

Additional  heave  for  6  ft. 
extra  fill 

HI 
H2 

0 
0 

Additional  pore  pressures  at 
failure  (in  ft.  of  water) 

P3 
P4 
P6 

27 

28.3 
7.9 
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TABLE  IV-14.   R.H.  BASSETT  -  PREDICTIONS 


QUANTITY 

PREDICTED 

VALUE 

Additional  height  of  fill  to 
cause  stability  failure 

13   ft. 

Additional  horizontal  move- 
ments for  6  ft.  extra  fill 

SI-3   El.  -30 

El.  -70 
SI-4   El.  -30 

El.  -70 

0.6"  1 

2.0"      All 

2.5"  [  eastwards 

2.25"J 

Additional  settlement  of  SP-1 
for  6  ft  extra  fill 

4.75" 

Additional  pore  pressures  (in 
ft.  of  water)  for  6  ft.  extra 
fill 

P3 
P4 
P6 

6 
8 
4 

Additional  heave  for  6  ft. 
extra  fill 

HI 

H2 

2" 
4" 

Additional  pore  pressures  at 
failure  (in  ft.  of  water) 

P3 
P4 
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FIGURE  4-31.   VERTICAL  DISPLACEMENTS  LOAD  INCREMENT  II 
Units  -  model  units  lxlO_1mm  =  1.8cm  (0.71  inch)  prototype 
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FIGURE  4-32.   HORIZONTAL  DISPLACEMENTS  LOAD  INCREMENT  II 
Units  -  model  units  lxl0_1mm  =  1.8cm  (0.71  inch)  prototype 
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V.   PRESENTATION  OF  FIELD  MEASUREMENTS 

5.1  Failure 

Early  in  the  morning  of  20  September  1974,  a  failure  of 
extraordinary  proportions  occurred.   Within  minutes,  a  simul- 
taneous failure  to  both  sides  of  the  embankment  caused  the 
crest  to  drop  about  30  feet  and  the  sides  to  heave  as  much 
as  14  feet.   An  aerial  view  of  the  slide  area  is  shown  in 
Figure  V-l.   Unfortunately,  no  one  witnessed  the  failure. 
No  surface  cracking  was  noticed  the  previous  day,  nor  was  a 
clear  indication  of  impending  failure  obtained  from  the  field 
instrumentation.   The  test  section  failed  at  an  elevation  of 
56.5  feet,  after  placing  18.7  feet  of  additional  fill  on  the 
existing  embankment. 

Additions  of  fill  to  the  test  embankment  had  proceeded 
much  according  to  plan.   The  first  three  construction  days 
were  used  to  steepen  the  east  slope  to  approximately  40°  by 
end-dumping  fill  from  the  embankment  crest.   Fill  was  moved 
with  trucks  and  the  embankment  crest  was  kept  level  with  a 
rubber-tired  loader.   Figure  V-2  shows  the  progress  of  fill- 
ing which  proceeded  at  an  average  rate  of  1.5  feet  per  day, 
excluding  weekends.   The  average  density  of  place  fill  was 
122  lb/ft   as  obtained  from  field  density  tests  during  fill- 
ing (see  Table  V-l) . 

5.2  Field  Measurements 
5.2.1   Pore  Pressures 

Figures  V-3  and  V-4  present  some  of  the  pore  pressure 
data  obtained  during  the  field  test.   The  computation  of  excess 
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pore  pressure  is  based  on  a  water  table  elevation  of  +  2.5  ft 
and  an  artesian  pressure  in  the  till  of  5  feet  of  water. 

5.2.2  Settlements  and  Heave 

Figure  V-5  is  a  plot  of  settlement  vs.  embankment 
elevation  for  the  four  settlement  plates  located  at  Station 
26  3.   Since  the  surveys  were  conducted  around  the  middle  of 
the  working  day,  the  embankment  elevation  selected  for  the 
settlement  and  heave  plots  was  the  average  elevation  during 
the  day  of  the  measurement. 

Figure  V-6  shows  the  vertical  movement  of  heave  points  H-l 
and  H-2  at  Station  263.   The  offset  for  H-l  was  130  feet  east 
and  for  H-2  183  feet  east. 

5.2.3  Horizontal  Movements 

Figures  V-7a  and  b  show  the  measured  horizontal  dis- 
placement at  inclinometers  SI-3  and  SI-4.   The  8  September  74 
plot  corresponds  to  movements  measured  two  days  after  an  addi- 
tional 6.6  ft  of  fill  had  been  in  place  for  two  days.   The 
arrows  in  the  figures  indicate  the  required  predictions. 

5.2.4  Summary  of  Measurements 

Table  V-2  contains  the  measurement  of  each  value  that 
the  predictors  were  requested  to  predict.   Also  shown  are  the 
values  predicted  by  each  team.   Two  values  are  given  for  the 
measured  pore  pressures  at  failure — the  values  obtained  the 
evening  preceeding  failure  and  the  values  obtained  one  to  two 
hours  after  failure. 
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5.3   Failure  Description 

A  surface  examination  immediately  after  the  failure  indi- 
cated cracks  and  movements  on  the  embankment  crest  and  sides 
from  Station  257  to  Station  267  +30.   Thus  failure  occurred 
for  a  length  of  1030  feet.   To  the  south  of  the  test  section, 
the  embankment  crest  was  at  its  original  40  ft  elevation  minus 
a  settlement  of  the  order  of  two  to  three  feet.   To  the  north, 
however,  the  40  ft  high  embankment  continued  only  to  Station 
266  after  which  the  embankment  elevation  decreased  to  +25  ft. 

Failure  occurred  to  both  sides  as  shown  in  Figure  V-8.   A 
single  frame  of  an  unsuccessful  time-lapse  film  made  during 
the  failure  indicates  failure  occurring  to  both  sides  simul- 
taneously and  movements  occurring  rapidly.   Also  shown  in 
Figure  V-8  are  the  depths  at  which  the  inclinometers  ruptured 
as  indicated  by  probing  after  the  failure  and  the  zones  of 
maximum  horizontal  strain  up  to  the  day  before  failure. 

At  Station  263,  the  top  of  the  embankment  dropped  a  dis- 
tance of  about  30  feet.   A  heave  of  14  feet  was  measured  on 
the  east  side  and  10  feet  on  the  west  side.   The  movements  to 
the  east  side  of  the  embankment  were  very  uniform  as  shown  in 
Figure  V-9.   A  very  large  crack  (14  ft  wide  at  the  top  and  8 
ft  deep)  formed  on  the  east  side  as  shown  in  the  photographs 
of  Figures  V-10  and  V-ll  and  schematically  in  Figure  V-9. 
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TABLE  V-l.   FIELD  DENSITY  TESTS 


TEST 
NO. 

DATE 

SAMPLE 
ELEV. 

DRY 
DENSITY 
lb/ft3 

WATER 
CONTENT,  % 

WET 
DENSITY 
lb/ft3 

1 

5  Sept  74 

44' 

118.0 

5.79 

124.5 

2 

5  Sept  74 

44' 

115.5 

4.92 

121.5 

3 

5  Sept  74 

44' 

102.5 

3.63 

108.5 

4 

6  Sept  74 

43.4* 

128.0 

7.05 

137.0 

5 

6  Sept  74 

43' 

127.5 

7.5 

137.5 

6 

6  Sept  74 

43' 

122.5 

3.95 

127.0 

7 

13  Sept  74 

48.7' 

102.0 

4.54 

109.0 

8 

13  Sept  74 

48.7* 

108.0 

2.57 

111.0 
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EXCESS  PORE  PRESSURE  VS  EMBANKMENT  ELEVATION 
STATION  263  TEST  SECTION 


•113- 


70 


60 


50 


40 


30 


20 


10 


P-6 
P-8 
P-9 
P-10 
P-11 

1  hr.  After  Failure 
7  hrs.  After  Failure 
—  Instrument 
Inoperative 


45  50  55 

EMBANKMENT  ELEVATION,  FT 


60 
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FIGURE  5-10. 
STATION  263  TEST  SECTION  AFTER  FAILURE  (East  Side) 


FIGURE  5-11.   LARGE  CRACK  140' 
EAST  OF  EMBANKMENT  CENTERLINE  AT  STA.  263 
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VI.   DISCUSSION  OF  RESULTS  AND  PREDICTIONS 
The  following  comments  were  made  by  the  predictors  immed- 
iately after  Dr.  Marr's  presentation  of  the  actual  field  re- 
sults.  During  this  period,  the  predictors  were  asked  to 
describe  their  immediate  thoughts  regarding  any  possible 
changes   or   alterations  they  would  make  in  their  prediction 
methods  based  on  the  actual  field  results.  These  presentations  were 
limited  to  approximately  three  to  five  minutes  each  and  the 
predictors  went  in  order  of  their  original  deliveries.   The 
first  predictor  was  Dr.  Harry  Poulos. 

Dr.  Poulos  would  have  used  the  same  analyses  with  one  ex- 
ception:  He  would  not  have  trusted  the  finite  element  stabil- 
ity analysis  as  much  as  he  had.   He  had  a  crude  prediction 
based  on  a  bearing  capacity  calculation  of  22  ft  which  was 
significantly  better  than  his  2  7  ft  submitted  prediction.   He 
felt  there  was  not  adequate  data  on  the  undrained  modulus  es- 
pecially stress-strain  time  properties.   He  would  have  liked 
to  have  seen  subsurface  earth  movement  measurements  in  order 
to  check  the  validity  of  the  overall  results  and  also  the  mech- 
anisms.  It  also  would  have  been  helpful  to  have  pore  pressure 
readings  in  the  till.   A  slightly  more  detailed  instrumentation 
set-up  would  have  avoided  the  situation  where  two  people  arrive 
at  the  same  answer,  which  is  the  correct  answer,  by  using  dif- 
ferent assumptions  and  different  methods,  and  having  no  way  of 
checking  out  which  one  is, in  fact, correct . 

The  next  predictor  to  speak  was  Bill  Moore  of  Dames  and 
Moore.   Mr.  Moore  was  pleased  to  see  some  attention  given  to 
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the  creep  phenomenon  which  has  bothered  him  for  quite  some 
time.   He  thinks  that  it  has  more  to  do  with  soil  behavior 
than  we  recognize  at  this  time.   He  raised  the  question  that 
if  the  fill  had  been  stopped  at  some  lower  elevation  and  then 
left  for  several  months,  would  it  in  fact,  have  failed.   Mr. 
Moore  felt  the  entire  exercise  was  most  constructive  especi- 
ally because  it  emphasized  the  requirement  for  good  data  and 
the  need  to  understand  soil  behavior.    Mr.   Moore  would 
like  to  see  laboratory  data  that  is  more  reflective  of  the 
actual  site  conditions.   He  questioned  some  of  the  field  pro- 
cedures that  were  not  well-documented.   For  example:   the 
settlement  of  SP1  and  whether  it  was  affected  by  the  removal 
of  peat  and  backfilling  which  apparently  was  not  well-docu- 
mented. 

The  next  predictor  to  speak  was  Professor  Charles  C.  Ladd. 
Professor  Ladd  said  that  he  would,  in  the  future,  stick  to 
heights  of  fill  and  deformations  and  avoid  comments  on  pore 
pressure.   Professor  Ladd  still  could  not  understand  why  the 
embankment  did  not  fail  at  a  higher  elevation. 

The  next  predictor,  Doctor  R.  H.  Bassett,  supported  Mr. 
Moore's  comment  regarding  the  strong  possibility  of  creep 
failure  had  the  embankment  -been  stopped  at  a  lower  elevation 
and  permitted  to  stand.   He  would  have  liked  to  have  had  in 
situ  block  samples  of  the  material  to  test  in  his  centrifuge. 
Boston  Blue  Clay,  with  its  sensitivity  of  3  to  5 ,  has  a  rather 
stiff  stress-strain  response  and  it  builds  up  rather  abruptly. 
When  the  clay  begins  to  fail,  we  find  that  finite  element 
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programs  can  not  deal  effectively  with  it.   Dr.  Bassett  is 
not  sure  in  fact,  whether  his  centrifuge  can  deal  effectively 
with  the  sensitive  clay.   He  is  fairly  sure  that  if  he  had 
built  the  embankment  to  a  height  of  15  ft  and  left  it  for  a 
period  of  nine  months,  it  would  have  failed  some  time  in  the 
period  of  0  to  9  months.   He  also  agreed  that  there  is  some- 
thing peculiar  regarding  creep,  especially  for  sensitive  clay, 
that  requires  more  serious  and  more  extensive  investigations. 

Dr.  Za-Chieh  Moh  of  the  Asian  Institute  of  Technology  would, 
if  he  had  to  do  it  over  again,  stick  to  a  simple  slip-circle 
analysis  but  put  more  emphasis  on  the  proper  soil  parameters 
and  the  proper  unit  weights  of  materials.   He  blames  MIT,  and 
particularly  Allen  Marr,  for  not  furnishing  adequate  data  for 
the  problem  especially  the  field  densities;  the  AIT  team  was 
lead  to  believe  that  it  was  much  lowor  than  122  lbs/cu  ft  and 
strength  data  was  also  inadequate.   He  stresses  very  much  that 
it  is  not  how  sophisticated  your  analyses  are  but  that  you  have 
accurate,  reliable  strength  data.   He  expresses  doubt  that  on 
a  real  practical  engineering  job,  he  would  have  enough  data 
to  justify  using  the  Finite  Element  Method  in  his  analysis. 
For  predicting  movements,  however,  he  feels  that  the  Finite 
Element  Method  is  the  only  method  to  use.   They  had  asked  for 
extension  test  data  and  for  creep  data  and  could  have  probably 
made  a  better  estimate  of  movements  had  they  had  this  data. 

The  next  speaker  was  Mr.  Verne  McGuffey  of  the  New  York 
State  Highway  Department.   Mr.  McGuffey  felt  that  with  enough 
study  of  the  data,  Professor  Ladd's  normalized  properties  con- 
cept and  the  limiting  equilibrium  analysis,  you  can  come  fairly 
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close  to  a  good  estimate  of  stability  for  this  problem.   He 
was  particularly  impressed  with  some  of  the  finite  element 
work  that  predicted  lateral  deformations.   This  work  could 
have  been  better  had  they  had  better  data.   He  also  felt  that 
the  unit  weight  of  material  especially  the  fill  in  place  per- 
haps was  a  little  low.   In  fact,  that  after  compaction,  it 
could  be  as  high  as  12  8  or  129  lbs/cu  ft.   For  the  settle- 
ment of  SP1  at  6  ft  of  fill,  they  used  an  extension  of  the 
settlement  that  had  taken  place  from  the  prior  6   ft    of 
filling.   After  studying  the  data  available,  they  concluded 
that  there  was  some  organic  material  left  in  place  below  the 
embankment  and  it  would  have  consolidated  during  the  period 
from  1969  to  1974.   Therefore,  under  the  new  loading,  it  would 
be  in  secondary  compression  and  would  behave  differently  from 
the  earlier  settlement  rate.   They  estimated  on  the  basis  that 
the  peat  was  still  in  virgin  compression,  which  was  in  fact, 
incorrect  and  they  plan  to  change  that  when  they  return  to 
their  offices. 

The  next  speaker  was  Professor  Duncan  from  the  University 
of  California  at  Berkeley.   His  comments  concerned  the  fact 
that  both  Professor  Ladd  and  he  were  using  Bishop's  Modified 
Method  and  essentially  the  same  way  of  evaluating  the  shear 
strength.   If  so,  then  why  did  Professor  Duncan  calculate  12 
feet  when  Professor  Ladd  calculated  23?   This  difference  must 
result  from  differences  in  the  details  of  the  calculations. 
Professor  Duncan  used  an  average  value  of  undrained  strength 
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whereas  Professor  Ladd  gave  much  more  weight  to  the  strength 
of  the  material  below  the  central  part  of  the  embankment. 
The  use  of  an  average  strength  had  seemed  like  small  detail 
when  he  was  making  the  prediction  but  in  fact,  on  second 
thought,  it  appeared  to  him  that  it  would  be  more  appropriate 
to  give  greater  weight  to  the  behavior  of  the  clay  beneath 
the  center  of  the  embankment.   Professor  Duncan's  second  point 
was  that,  if  you  take  the  extreme  heights  of  fill  at  failure 
predicted,  the  actual  variation  is  only  plus  or  minus  20% 
from  the  actual  failure  elevation.    He  suggested  that  this 
was  a  more  realistic  way  of  looking  at  the  numbers,  rather 
than  comparing  the  additional  height  of  fill.   He  expressed 
the  opinion  that  our  predictive  ability  for  stability  is  in 
fact  quite  good,  and  he  congratulated  MIT  on  both  their  test- 
ing and  their  ability  to  perform  analyses.   He  again  emphasized 
the  importance  of  creep,  especially  relative  to  movements. 
Had  he  waited  a  length  of  time  with  the  six  feet  of  fill  in 
place,  the  movements  would  have  been  considerably  larger.   We 
don't  know  how  to  scale  creep  effects  from  the  lab  to  the  field 
and  we  don't  know  how  to  subtract  the  creep  effects  out  of 
measured  movements.   Creep  is  an  important  factor,  not  atfi'ly**  - 
in  stability  calculations  but  in  analyzing  field  movements  and 
in  interpreting  lab  test  results  and  projecting  them  to  the 
field. 

The  next  speak,  Mr.  Robert  C.  Kirby,  supported  Prof.  Dun- 
can's observation  that  the  differences  in  predictions  were  not 
large  when  one  considers   the  total  height  of  fill  to  cause 
failure  rather  than  the  incremental  height  of  fill  to  cause 


126- 


failure.   He  believes  that  the  most  important  result  of  the 
failure  is  a  back-calculated  undrained  shear  strength  for 
Boston  Blue  Clay.   Mr.  Kirby ' s  analyses  (and  those  of  other 
predictors)  showed  one  could  predict  a  wide  range  of  failure 
heights  depending  on  what  shear  strength  profile  is  used. 

There  is  a  difference  in  philosophy  between  analysis  and 
design.   We  were  not  asked  to  design  an  embankment  on  soft 
clay  but  rather  to  predict  the  failure  of  an  embankment  on 
clay.   For  design  of  this  embankment,  one  would  probably  de- 
termine shear  strength  from  unconfined  compression  tests  and 
estimate  the  design  elevation  of  the  embankment  to  be  equi- 
valent to  a  loading  of  about  six  times  the  undrained  shear 
strength  using  a  factor  of  safety  of  1.1.   If  the  embankment 
had  to  be  higher  than  this,  one  could  make  field  vane  tests 
in  order  to  better  determine  shear  strength.   However,  the 
design  elevation  of  the  embankment  would  be  determined  using 
a  higher  factor  of  safety.   If  one  were  designing  for  sta- 
bility based  on  consolidated,  undrained,  triaxial,  compression 
tests,  an  even  higher  factor  of  safety  would  be  used  when 
arriving  at  the  design  fill  elevation. 

When  hindcasting  measured  performance  using  the  Finite 
Element  Method,  one  can  "predict"  almost  anything  depending 
on  the  parameters  which  are  input  to  the  model.   There  are  at 
least  eight  or  ten  parameters  which  can  be  adjusted  and  the 
answers  are  controlled  by  the  mix  of  parameters  one  chooses. 
Mr.  Kirby  feels  that  some  of  the  judgments  made  relative  to 
selection  of  parameters  for  finite  element  solutions  were  not 
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reasonable  because  they  were  not  supported  by  the  available 
data. 

The  next  predictor  that  spoke  was  Professor  Gerald  Raymond 
from  Canada.   He  said  that  the  vane  strengths  were  lower  on 
the  west  side  than  on  the  east  side  and  so   he  did   feel 
it  could  fail  both  ways.   Professor  Raymond  would  have  placed 
a  lot  more  emphasis  on  vane  tests.   He  feels  that  lateral 
movements  are  important  and  should  be  measured  and  understood 
and  you  have  to  go  to  a  finite  element  analysis  to  deal  with 
these  horizontal  movements.   He  is  still  modifying  his  model 
which  is  similar  to  the  Cambridge  model.   Professor  Raymond 
feels  that  if  you  are  going  to  put  numbers  into  a  computer 
program,  you  may  as  well  use  a  complex  computer  program. 

The  next  predictor  to  speak  was  Doctor  Peter  Wroth  from 
England  who  indicated  that  some  of  his  colleagues  thought 
he  was  mad   to  get  involved  in  such  a  symposium.   He  compli- 
mented MIT  on  the  conduct  of  the  entire  effort  especially  the 
actual  causing  of  the  failure  within  the  specified  time  limit. 
Dr.  Wroth  said  that  from  the  engineering  standpoint,  there  are 
two  processes  involved.   One  is  prediction  and  design  and  the 
other  is  analysis  and  understanding:   they  are  distinctly  dif- 
ferent processes.   Regarding  the  first  of  these,  they  had  used 
elastic-plastic  models  of  soil  behavior  and  were  particularly 
interested  in  the  predictions  of  excess  pore  pressures.   He  was 
very  satisfied  with  his  finite  element  predictions.   He  felt  it 
would  be  interesting  to  ask  the  predictors  to  resubmit  predic- 
tions of  movements  and  pore  pressures  for  the  actual  height  of 
the  embankment  at  failure  rather  than  the  (different)  heights 
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that  each  had  predicted  would  cause  failure.   Regarding  the 
second  process  of  analysis  and  understanding,  he  considered 
that  we  should  now  see  if  any  change  of  behavior  of  the  soil 
from  elastic  to  plastic  had  occurred  prior  to  failure.   This 
could  be  done  by  a  careful  evaluation  of  some  of  the  field 
measurements . 

Next,  Prof.  Lambe  asked  the  predictors  to  comment  on  whether 
they  were  surprised  at  the  extent  of  the  failure,  in  other  words, 
its  propagation  both  north  and  south  from  the  loaded  area. 

Dr.  Wroth  was  surprised!   He  said  that  most  methods  of 
analysis  were  two-dimensional  and  did  not  consider  the  effect 
of  the  length  of  propagation.   If  asked  whether  the  sliding 
mass  would  propagate  beyond  the  loaded  area,  he  would  have 
said  yes,  but  certainly  nothing  like  the  extent  to  which  it 
had  happened. 

Prof.  Duncan  felt  that  the  extent  of  the  failure  was  one 
of  the  big  mysteries  of  the  entire  symposium.   He  is  not  sure 
why  it  failed  in  such  a  fashion. 

Dr.  Bassett  was  not  surprised  at  the  extent  of  the  failure 
and  again,  postulated  that  the  very  high  pore  pressures  gener- 
ated within  the  failure  zone  could  probably  be  transmitted 
quite  rapidly  in  a  horizontal  direction  and  cause  additional 
failures. 

Dr.  Wroth  felt  that  some  work  should  be  done  to  try  to  esti- 
mate the  amount  of  energy  that  was  dissipated  as  the  failure 
propagated  away  from  the  loaded  zone.   Is  this  eneryg  trans- 
mitted by  means  of  the  high  excess  pore  pressures,  or  is  it 
transmitted  by  means  of  the  intermediate  effective  stress,  a2? 
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Professor  Raymond  was  concerned  that  the  embankment  would 
fail  in  both  directions  and  felt  vindicated  by  the  results. 

Mr.  McGuffey  described  the  recent  failure  of  about  1,000 
ft  of  interstate  embankment  that  propagated  an  additional 
2000  ft.   One  of  the  safety  factors  in  the  failure  zone  was 
about  .95  where  is  was  overstresses  and  calculated  1.4  in  the 
other  zones  prior  to  failure.   He  suggested  that  the  high- 
pore  pressure  is  generated  at  the  failure  zone   and  can,  in 
fact,  travel  quickly,  rapidly,  and  at  great  distances  in  a 
lateral  direction.   Prof.  Raymond  feels  that  the  sensitivity 
of  the  clay  had  something  to  do  with  the  propagation  of  the 
failure  zone  beyond  the  loaded  zone.   Dr.  Bassett  felt  that 
perhaps  the  clay  is  laminated  and  then  can  transmit  these  pore 
pressures  laterally  very  quickly.   Professor  Duncan  said  that 
he  sees  very  little  use  for  predicting  pore  pressures  for  this 
type  of  problem — that  movements  can  be  predicted  and  stability 
analyzed  without  predicting  pore  pressures. 

Professor  Ladd  said  that  an  effective  stress  analysis  at 
failure  is  absolutely  worthless. 

Dr.  Wroth  emphasized  that  his  finite  element  analysis  is 
an  effective  stress  analysis.   Since  it  is  effective  stresses 
that  control  the  behavior  of  the  soil,  it  is  those  stresses 
that  we  should  be  studying.   Dr  Wroth  said  that  as  long  as  you 
are  studying  pre-failure  conditions,  it  was  proper  and  appro- 
priate to  use  effective  stresses,  but  that  conditions  at  and 
post  failure  might  be  better  predicted  in  terms  of  total  stresses 
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because  of  the  difficulties  of  predicting  post-failure  pore 
pressures. 

Prof.  Moh  also  agrees  that  the  lateral  extent  of  a  failure 
is  related  to  pore  pressure  transmission,   Mr.  Kirby  said  that, 
based  on  his  experience,  limiting  equilibrium  analyses  using  measured  pore 
pressures  and  effective  stress  parameters  do  not  give  satis- 
factory predictions.   Mr.  Moore  also  agrees  with  Prof.  Duncan 
in  that  pore  pressures  are  rather  useless  to  understand  the 
failure  conditions. 

Dr.  Bassett  also  felt  that  his  increase  in  pore  pressure 
did  propagate  and  did  in  fact,  contribute  to  the  failure  on 
either  side  of  the  loaded  section.   Professor  Raymond  agrees 
with  that  statement  although  he  says  creep  is  also  very  im- 
portant in  this  regard.   Mr.  McGuffey  said  that  he  does  believe 
in  creep  and  that  they  have  had  progressive  failures  in  New 
York.   Failures  have  taken  place  approximately  a  month  after 
the  end  of  filling.   Calculated  factors  of  safety  were  about 
1.3.   Mr.  McGuffey  says  that  they  run  special  tests  in  New  York 
on  strains  on  the  order  of  20%  and  can  modify  normal  stress-strain 
relationships  when  they  find  areas  in  an  embankment  analysis 
that  become  overstressed.   These  areas1  strength  is  then  reduced 
before  performing  the  stability  analysis  and  thus  a  more  real- 
istic value  of  safety  factor  is  obtained  by  using  these  reduced 
strengths  in  areas  that  have  been  overstressed. 

At  this  point,  the  discussion  focused  on  the  differences 
between  and  importance  of  creep  and  progressive  failure.   Dr. 
Wroth  said  that  his  understanding  of  creep  is  deformation  that 
occurs  under  a  constant  effective  stress  whereas  for  progressive 
failure,  the  effective  stresses  are  changing. 
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Professor  Raymond  said  that  undrained  creep  would  in 
fact  generate  pore  pressures;  therefore,  the  effective  stress 
would  not  be  constant.   Mr.  Moore  said  that  the  issue  is  a 
very  important  one  as  time  dependent  creep  may  occur  at  very 
low  stresses.   So,  Mr.  Moore  says  we  are  novv  again  worrying 
about  stress  levels  that  could  cause  slow  movements  that  would 
then  continue  on  to  failure  which  is  sort  of  the  same  thing 
that  we  worried  about  30  years  ago  under  the  name  of  yield 
point  stress. 

Professor  Lambe  then  asked  the  panel  how  they  would  con- 
trol construction  if  they  wanted  to  build  this  embankment  as 
quickly  as  possible  and  as  high  as  possible  while  keeping  the 
factor  of  safety  very  slightly  over  1.  Professor  Ladd  would 
use  inclinometer  measurements  to  control  construction  because 
he  feels  that  lateral  deformations  are  the  best  indicator  of 
impending  failure. 

Professor  Duncan  says  that  if  you  ran  a  total  stress  analy- 
sis, then  you  must  use  movements  to  control  construction  and 
he  feels  that  movements  are  a  very  difficult  way  to  do  this. 
The  reliability  of  movement  predictions  is  far  less  than  the 
reliability  of  predictions  of  stability.   Mr.  Moore  said  that 
he  feels  movement  measurements  are  very  important  and  they  do 
this  frequently.   Especially,  he  looks  at  the  relationship  of 
volume  change  due  to  vertical  movement  and  horizontal  move- 
ment.  They  use  a  combination  of  vertical  movements  and  hori- 
zontal movements  to  note  whether  they  are  getting  consolida- 
tion and  increase  in  strength  in  addition  to  other  phenomena. 
Mr.  Moore  also  said  it  is  important  to  evaluate  the  time  rate 
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of  movement.   Mr.  Kirby  feels  that  controlling  construction 
with  field  measurements  is  very  difficult;  there  is  no  one 
field  measurement  which  conclusively  indicates  imminent 
failure.   Measurements  of  settlement  and  horizontal  movement 
are  useful  indicators.   In  addition,  he  would  recommend  pore 
pressure  measurements.   A  plot  of  total  head  vs.  fill  ele- 
vation should  be  made  for  each  piezometer.   A  change  in  slope 
of  the  data  from  one  straight  line  segment  to  another  indi- 
cates the  soil  at  that  location  has  yielded.   We  should  not 
be  concerned  when  the  soil  beneath  the  center  of  the  embank- 
ment yields  locally  but  we  are  concerned  when  the  yielded 
zone  approaches  the  toe  of  the  embankment.   Consequently,  pie- 
zometers should  be  installed  near  the  toe  of  the  embankment. 
Dr.  Moh  would  definitely  use  piezometers  to  monitor  pore  pres- 
sures and  an  effective  stress  analysis  and  therefore  to  keep 
calculating  stability  as  the  embankment  grew.   Dr.  Moh  says 
that  if  you  measure  lateral  movements,  by  the  time  you  got 
significant  ones,  it  would  probably  be  too  late. 

Professor  Lambe  asked  the  panel  whether  any  of  them  would 
have  or  could  have  used  decision  theory  or  probability  in 
coming  up  with  their  prediction.   Dr.  Wroth  was  the  first  to 
respond.   Before  we  could  adopt  such  theories,  we  must  have 
some  parameters  that  we  have  confidence  in,  or  a  certain  degree 
of  confidence.   Dr.  Wroth  was  grateful  to  MIT  for  exposing  him 
to  the  need  to  make  use  of  probability  and  decision  theory  and 
he  would  make  an  effort  to  consider  it  in  detail. 
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Mr.  McGuffey  felt  that  you  could  not  apply  theories  of 
probability  and  risk  analysis  in  a  case  where  life  is  in 
danger  that  in  fact,  there  you  must  have  100%  reliability. 

Mr.  Moore  felt  that  in  a  typical  project,  you  rarely  have 
enough  data  to  apply  statistical  probability  theory.   Mr, 
Moore  also  felt  that  in  engineering,  there  is  no  such  thing 
as  100%  safety. 

Professor  Ladd  had  a  comment  regarding  probability  in 
that  there  is  some  bias  introduced  because  of  the  nature  of 
the  test  data.   One  can  run  several  kinds  of  strength  tests  and 
get  very  different  answers.   We  then  have  to  be  able  to  resolve 
the  fact  that  there  is  uncertainty  in  the  measured  strengths . 
So,  relating  this  further  to  probability  of  failure  becomes 
very  difficult. 

Dr.  Bassett  had  a  comment  that  stems  from  a  job  he  did  in 
the  United  Kingdom  where  he  installed  500  soil  anchors  and 
tested  them  all  to  1.5  times  working  stress  and  failed  some 
3  9  of  these  anchors.   His  standard  deviation  band  was  much 
wider  than  he  had  predicted.   It  very  strongly  suggested  to  him 
that  on  a  job  where  there  are  a  large  number  of  repetitive 
operations  of  a  similar  nature,  that  we  cannot  very  easily 
decide  what  is  or  what  is  not  reasonable  to  expect  for  a  stan- 
dard deviation.   He  also  added  that  use  for  probability  should 
probably  be  restricted  to  jobs  that  are  much  like  mechanical 
engineering  where  we  have  this  repetitive  action  that  we  can 
apply  the  theory  of  probability  to. 
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Prof.  Raymond  was  not  very  enthusiastic  about  statistics 
and  probability  again  claiming  that  the  primary  reason  is  the 
lack  of  sufficient  data  to  call  it  a  representative  sample. 

After  the  discussion  on  decision  theory,   Prof.  M.  M. 
Baligh  of  MIT  discussed  the  three-dimensional  aspects  of  the 
failure.  In  the  initial  formulation  of  the  problem  itself,  Dr. 
Marr  had  considered  the  length  of  the  loaded  section  in  order 
to  achieve  a  plane  strain  situation.   The  extent  of  the  failure 
beyond  the  limits  of  the  loaded  area  indeed  surprised  Prof. 
Baligh.   Prof.  Baligh  reviewed  the  literature  and  found  very 
few  three-dimensional  slope  stability  analyses  that  responded 
to  the  needs  of  this  particular  problem.   Prof.  Baligh  considered 
several  geometric  shapes,  simple  cylindrical  slip  surface  with 
straight  ends  and  then  with  elipsoid  and  conical  shaped  ends 
trying  to  find  the  variations  with  safety  factor  versus  some 
of  his  geometric  parameters.   Prof.  Baligh  concluded  that  the 
end  effects  increase  the  factor  of  safety  obtained  by  means  of 
a  conventional  2-D  solution  by  about  19%  for  the  1-95  test  fill. 
In  addition,  he  said  that  even  though  the  end  effects  on  the 
factor  of  safety  can  be  predicted  with  reasonable  accuracy,  the 
length  of  failure  is  very  sensitive  to  the  parameters  defining 
the  problem  and  is  thus  difficult  to  predict. 

After  Prof.  Baligh 's  presentation,  there  were  questions 
for  the  panel  from  the  floor.   These  questions  are  briefly  sum- 
marized below. 

Prof.  Selig  asked  that  if  they  were  building  this  embank- 
ment from  scratch,  starting  from  ground  zero  and  building  up, 
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would  they  have  used  the  vane  values  and  had  confidence  in 
them?   Dr.  Poulos  said  he  would  not  use  the  vane  shear  values 
without  trying  to  correlate  them  with  some  type  of  lab  data. 
Professor  Duncan  said  that  he  was  very  surprised  that  the  UU 
tests  gave  strengths  that  were  only  40%  of  the  strength  de- 
termined using  CU  tests  for  Boston  Blue  Clay.   Professor  Dun- 
can said  on  a  design  job,  he  might  use  the  vane  shear  strength. 
Mr.  Kirby  responded  that  if  you  had  no  prior  experience  in  the 
area  and  only  vane  shear  strength  values,  you  would  have  to  be 
satisfied  using  a  rather  high  safety  factor,  perhaps  as  high 
as  1.6.   Prof.  Raymond  said  that  in  Canada  they  put  a  lot  of 
faith  in  vane  strength  values,  design  many  embankments  which 
in  fact  stand  up  if,  of  course,  the  vane  test  is  run  properly. 
He  added  that  many  engineers  prefer  a  field  test  over  a  lab 
test.   Dr.  Wroth  thought  that  sampling  disturbance  was  always 
a  problem  and  that  in  particular,  the  unconfined  compression 
test  was  a  complete  waste  of  time  and  effort. 

Mr.  McGuffey  basically  doesn't  use  the  U  and  UU  tests.   He 
feels  that  you  can  get  as  much  information  and  more  reliable 
and  consistent  information  with  the  pocket  penetrometer.   Mr. 
Kirby  was  asked  a  question  about  safety  factor  chosen  as  a  func- 
tion of  the  type  of  strength  data.   He  pointed  out  that  most 
of  our  experinece  is  based  on  U  and  UU  tests  and  that  the  fac- 
tor of  safety  adopted  for  design  is  based  on  case  histories 
where  U  and  UU  tests  were  used.   More  recently,  we  are  attempting 
to  use  some  of  the  newer  tests  that  give  different  (and  usually 
higher)  shear  strengths.   We  don't  have  a  long  history  of  exper- 
ience with  these  test  results  and  must  reconsider  the  safety 
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factor  that  we  use.   In  his  view,  a  higher  safety  factor  is 
required  for  design  when  using  the  more  sophisticated  labora- 
tory tests.   Prof.  Ladd  reminded  those  present  that  throughout 
the  world,  we  don't  always  get  a  lower  strength  value  with  a 
UU  test.   The  Boston  Blue  Clay  samples  we  used  were  taken  from 
very  large  depths  and  therefore  we  had  much  lower  strengths 
due  to  excessive  disturbance. 

The   question  was   asked    "if   they   had   to   design   another   em- 
bankment  adjacent   to   this   one,    what    safety   factor   and  method 
would   they   use?"      Dr.    Bassett   responded   that   he  would   use    a 
Bishop's    analysis  with   vane    shear   strengths   and   no   correction 
factor.      Mr.    McGuffey  would   have   used   a   safety   factor   of    1.25 
with   the    same    strength   parameters.      Dr.    Wroth  would   use   vane 
shear    strengths   but   he  would   ask   the   client  what  the   conse- 
quences  of   failure  would  be   before   choosing   a   safety   factor. 
Prof.    Raymond  would   have   used   the   Bishop   analysis,    vane    shear 
strengths    and   a    safety    factor   of    1.15.      Professor   Duncan    said 
that   he  would   use   the   SHANSEP    stress   parameters   with   a   creep 
correction  and   a   safety   factor  of   1.15.      Dr.    Poulos  would  have 
used   a   Bishop   analysis,    SHANSEP    strength  parameters   with   a   cor- 
rection  for   anisotropy   and   a    safety   factor  of    1.3.      Mr.    Kirby 
would  run  field  vanes  along  the  alignment  of  the  new  embankment  and  if 
they  were  no  lower  than  the  previous  one,   he  would  not  do  any  analysis 
and  would  tell  the  client  he  could  build  to  elevation  40  with  no  trouble. 
A  question  was  asked  from  the  audience:      "Is  there  any  way  to  pre- 
dict the  increase  in  pore  pressure  outside  the  failure  zone?"     Mr.   Kirby 
said  that  he  doesn't  feel  that  pore  pressures  transmitted  laterally  caused 
the  additional  failure.     He  feels  that  the  ends  of  the  embankment  were  just 
dragged  down  by  shear  stresses  and  then  the  failure  propagated  farther 
and  farther  away  from  the  original  failure  area.     Prof.  Raymond  reminded 

people  that  the  actual  loaded  area  was  probably  too  short  to  be  considered 
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plane  strain.   There  was  a  very  strong  possibility  that  the 
initial  failure  would  have  been  experienced  at  a  lower  height 
than  the  18.5  feet  which  caused  the  failure  at  our  particular 
case.   Robert  Tierney  of  the  MDPW  asked  the  panel  "If  we  de- 
cide that  we  are  going  to  build  1-95  and  we  now  need  this 
road,  what  do  we  do  with  our  completely  remolded  foundation?" 
Prof.  Raymond  said  that  it  will  get  large  differential  settle- 
ments in  this  area  because  it  is  remolded  and  it  has  a  heavier 
load  on  top  of  it.   Mr.  McGuffey  says  that  he  agreed  that  it 
would  get  large  differential  settlements  and  therefore  he  would 
place  several  berms  and  put  another  surchage  load  on  the  area 
to  try  to  even  up  some  of  these  zones  and  build  a  highway.   Dr. 
Bassett  felt  that  he  wouldn't  dare  put  a  surcharge  load  on  this 
area  and  that  he  would  build  a  road  and  just  plan  to  resurface 
it  every  six  months  or  so  until  the  large  settlement  stopped. 
Dr.  Marr  felt  that  even  after  the  failure,  we  are  still  above 
the  elevation  of  the  future  road  and  that  if  it  were  not  built 
up  for  six  months  to  a  year,  we  would  actually  be  taking  material 
off  at  that  time  and  have  had  some  further  consolidation  and 
settlement  that  would  help  us. 

A  question  was  asked  by  an  engineer  from  TAMS ,  Ernest  Jonas, 
whether  the  strength  of  the  fill  was  considered  in  the  analyses, 
since  there  was  reasonable  height  of  fill  and  shearing  took 
place  in  the  fill.   Prof.  Raymond  said  that  he  did  not  consider 
the  strength  in  his  analyses.   Joe  Olsen  of  MIT  asked  the  ques- 
tion "If  we  had  this  embankment,  and  the  problem  was,  instead  of 
trying  to  fail  it  to  put  on  an  additional  height  of  fill  that 
was  to  last  for  a  minimum  of  2  5  years,  if  we  considered  creep 
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phenomena,  what  is  the  height  of  fill  that  the  panel  would 
agree  to  place  on  it  if  it  had  to  last  25  years?"   Mr  McGuffey 
said  13  feet   Dr.  Basset  said  13  feet  also.   Prof.  Raymond 
would  go  between  8  and  13  feet.   Dr.  Wroth  said  that  he  could 
pick  14  ft  out  of  the  air,  but  that  the  phenomenon  of  creep 
raised  very  complex  problems  that  should  be  looked  at  in  detail 
Prof.  Duncan  said  that  he  feels  that  the  creep  problem  does 
go  quickly  so  he  probably  wouldn't  add  more  than  about  5  feet 
of  fill  to  the  embankment.   Prof.  Ladd  would  backfigure  the 
strength  for  the  18.7  feet  and  then  use  a  safety  factor  of  1.3 
to  reduce  the  height.   Dr.  Moh  would  have  put  on  about  8  or  9 
feet.   Mr.  Kirby  said  he  would  not  know  how  to  answer  that  ques- 
tion.  Mr.  Kirby  is  not  convinced  that  we  have  a  creep  problem 
with  Boston   Blue  Clay.   He  can  easily  predict  the  amount  of 
consolidation  settlements  that  are  required  to  fit  the  problem; 
therefore,  in  his  opinion,  creep  is  not  a  major  problem  here. 

Professor  Duncan  and  Mr.  Kirby  got  into  a  discussion  about 
how  much  of  the  time-dependent  settlement  was  due  to  consolida- 
tion and  how  much  was  due  to  creep. 

Mr.  McGuffey  agrees  with  Professor  Duncan  in  that  we  can- 
not explain  all  of  the  settlements  by  using  consolidation  set- 
tlements.  The  discussion  was  cut  off  at  this  point  and  Prof. 
Lambe  delivered  a  short  summary  prior  to  the  closing  of  the 
symposium. 
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VII.   CONCLUDING  REMARKS 

Professor  Lambe  felt  that  the  symposium  turned  out 
very  well.   He  said  that  we  always  have  good  presentations 
when  we  get  able  people  to  focus  on  a  single  problem. 

We  started  out  with  the  premise  that:   predictions 
constitute  an  important,  a  very  essential  part  of  Civil 
Engineering.   Making  the  prediction  is  very  necessary,  but 
insufficient.   The  selection  of  parameters  is  very  important 
and  is  the  major  limitation  in  making  good  predictions. 
That  the  accuracy  of  the  prediction  is  directly  related  to 
the  sophistication  of  the  analysis  is  a  totally  erroneous 
statement.   We  should  leave  this  symposium  understanding 
that  fact. 

There  are  three  classes  of  predictions:   the  crude 
prediction,  the  engineering  prediction  and  the  state-of-the 
knowledge.   We  do  not  necessarily  get  a  more  accurate  pre- 
diction if  we  use  a  state-of-the-knowledge  prediction. 
There  are  two  reasons  for  making  the  state-of-the-knowledge 
prediction:   first,  it  gives  an  insight  into  the  fundamentals 
of  the  problem  and  second,  it  gives  a  good  basis  for  a  future 
evaluation.   It  is  very  difficult  to  use  the  crude  analysis 
and  then  make  an  evaluation  of  the  prediction. 

Prof.  Lambe  sees  the  finite  element  technique  as  being 
a  very  powerful  tool  and  playing  two  important  roles:  1)  to 
generate  parametric  studies  that  can  look  at  a  problem  in 
several  different  ways  and  2)  to  provide  a  very  powerful 
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analysis  tool  for  the  state-of-the-knowledge  type  analysis. 
He  feels  that  the  average  engineer  will  not  have  the  time 
or  money  to  use  a  finite  element  analysis  on  the  first  cut. 

Prof.  Lambe  was  impressed  with  the  centrifuge  test. 
However,  he  felt  its  main  use  will  be  identifying  mechanisms 
and  not  solving  practical  problems  by  making  numerical  pre- 
dictions. 

He  was  impressed  with  the  extrapolation  technique 
which  is  the  simplest  way  to  make  a  prediction.   We  can  extra- 
polate to  get  stresses,  movements,  pore  pressures,  but  we 
cannot  extrapolate  very  well  to  get  factor  of  safety.   Ex- 
perience is  very  important  but  probably  one's  own  experience 
has  the  most  significance.   Sometimes  in  engineering  it  is 
rather  difficult  to  read  of  the  experiences  of  others  and 
then  relate  to  them  accurately.   For  experience  to  be  use- 
ful, it  must  be  evaluated  experience.   If  you  have  done  the 
same  thing  over  and  over  again  for  twenty  years,  you  have 
twenty  years  *  experience ,  but  you  are  probably  not  much  bet- 
ter off  than  you  were  at  the  end  of  the  first  year. 

The  symposium  has  raised  a  number  of  important  ques- 
tions and  Prof.  Lambe  felt  that  as  we  think  about  these 
questions  hopefully  we  can  come  up  with  better  explanations 
than  those  offered  at  the  sumposium.   Why  the  embankment 
failed  over  such  a  large  area  was  not  answered  properly. 

Prof.  Lambe  concluded  the  symposium  with  thanks  to 
the  predictors  for  the  preparation  and  presentation  of  their 
predictions. 
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Richard  H.  Bassett 
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Dr.  Bassett  graduated  from  King's  College,  London  in  1958. 
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responsible  for  the  foundation  work  and  design  of  several  major 
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Asian  Institute  of  Technology 

Dr.  Brand  received  his  B.  Sc .  and  Ph.  D.  from  Leeds  Univer- 
sity.  He  is  the  Editor  of  the  journal,  Geotechnical  Engineering. 
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Dr.  Davis  was  educated  at  University  College  School  and  Uni- 
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Laboratory  in  the  UK  and  a  lectureship  at  University  College, 
he  joined  the  staff  on  the  University  of  Sydney  in  1952,  becom- 
ing a  Professor  of  Civil  Engineering  in  1968.   He  is  a  Fellow 
of  the  Institution  of  Engineers,  Australia  and  is  one  of  the 
founders  of  the  Australian  Geomechanics  Society.   He  is  active 
in  the  affairs  of  the  International  Society  for  Soil  Mechanics 
and  Foundations  having  been  a  Vice-President  of  that  Society. 
He  has  spent  periods  of  study  leave  at  MIT,  the  Building  Research 
Station  in  the  UK  and  Imperial  College,  London.   The  two  main 
areas  of  his  research  are:   the  development  of  a  unified  theory 
of  the  behavior  of  foundations  and  the  use  of  the  theories  of 
elasticity  and  plasticity  for  soil  and  rock  stability  problems - 

James  Michael  Duncan 
Department  of  Civil  Engineering 
Unviersity  of  California,  Berkeley 
Berkeley,  California 

Professor  Duncan  received  his  BSCE  in  19  59  and  his  MSCE  in 
1962  from  the  Georgia  Institute  of  Technology.   In  1965,  he  re- 
ceived his  Ph.  D.  from  the  University  of  California  at  Berkeley. 
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He  joined  the  staff  at  the  University  of  California,  Berkeley  in 
1964.   Professor  Duncan  is  currently  Chairman  of  the  Committee  on 
Embankment  Dams  and  Slopes  of  the  American  Society  of  Civil  Engi- 
neers.  In  1972,  he  was  the  recipient  of  the  ASCE  Collingwood 
Prize,  and  in  19  73  he  was  awarded  the  Walter  Huber  Research  Prize 
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modeling  of  stress-strain  behavior  of  soil,  and  analysis  of 
slope  stability  problems.   He  is  a  registered  Civil  Engineer 
in  California. 
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